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FOREWORD 


This volume contains the programming documentation for the 
NASTRAN updates and the MESHGEN program for the NASTRAN Hydroelastic 
System. Volume I contains the Theory and Results. Volume III 
contains the User Manual updates. 
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4.0 NASTRAN PROGRAM MANUAL UPDATES 


The operational steps, data descriptions, and program code for the new 
NASTRAN hydroelastic analysis syste.^i are described in this section. The format 
is similar to that of the NASTRAN Programmer's Manual for the module and sub-' 
routine descriptions. The overall f?ow of the system is described, followed 
by the descriptions of the individual modules and their subroutines. 

4.1 OVERALL PROGRAM FLOW 

From a programmer's viewpoint the NASTRAN hydroelastic normal modes anal- 
ysis is a straight forward extension of Rigid Format 3, the structural normal 
modes aialysis. The overall flow of the problem for a direct formulation, 
shown itt Figure 1, describes the key steps and data blocks used in the DMA? 
program. Because of the special characteristics of the fluid degrees of free- 
dom, the USET definitions in NASTRAN have been modified. These components 
describe the size and characterxstics of the matrices and vectors. 

The names of these new degree-of-freedom sets are defined below: 

u : All structure and fluid components 

g 

u^: Structure components 

Uyt Fluid components 
u^^: Free surface degrees of freedom 

u^: Reduced structure components 

u^: Reduced structure plus free surface displacements = solution set 

In the DMAP the names of matrices and vectors generally correspond to these 
sets. For instance, AXY is a matrix with rows defined by structure points (x) 
and columns corresponding to fluid points (y) . 

The steps described in Figure 1 are: 

Step 1: The normal steps used in Rigid Format 3 are used to build the stiff- 

ness (KGG) and mass (MGG) matrices for both structu.re and fluid. 

The matrix terms for the fluid points are actually pressure/ flow 
functions rather than the normal displacement/force functions. The 
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FIGURE 1. NASTRAN HYDROELASTIC DMAP FLOW (Direct Assembly) 
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FIGURE 1. NASTRAN HYDROELASTIC DMAP FLOW (Cont’d) 
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FIGURE 1. NASTRAN HYDROELASTIC DMAP FLOW (Cont'd) 
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Step 2 


Step 3 


Step 4 


USET data block, generated by module GP4, describes all degrees of 
freedom. 

; The FLBMG (Fluid/Boundary Matrix Generator) Module generates the area 
factor matrix, AF, and the additional stiffness due to gravity, DKGG. 
It also produces two new USET data blocks. USETF contains the fluid 
point flags and USETS describes only structure points. 

: The following DMAP partitioning operations are coded in the system; 

a) Partition g -»• x + y 
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b) Partition y -*■ fr + j 
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The SWITCH module is used to replace the USET data block with the 
new block which describes only structure points. The normal struc- 
ture processing modules will use this data block and operate only 
on the subsets of the u^ set. 
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Steps 5 & 6: For non-restart problems, the structure stiffness and mass 

matrices are constrained and reduced In the usual manner. For 
restart problems with no structure changes, these operations are by- 
passed by the use of a PARAM bulk data card. 

Step 7: For fluid prob'.ems with gravity effects, the incremental stiffness 

matrix, DKXX, is constrained and reduced to correspond to the reduced 
structure points. 

Step 8: Module GFSMA (General Flu Id /Structure Matrix Assembler) performs the 

operations to merge the structure and fluid matrices Into the solu- 
tion matrices (see Section 4.4 for the detailed steps). The output 
matrices correspond to the reduced structural displacements and the 
free surface displacements. 

Step 9: The NASTRAN solution modules are executed at this stage. For normal 

modes, the READ (Real Eigenvalue Analysis - Displacement method) 
module is used. However, the capability exists to restart In a 
dynamics analysis rigid format. 

Step 10: MPYAD Is used to recover fluid point pressures and free surface point 
displacements: 


The structure and free surface displacements are separated, d -*■ a -f fr. 


m -> 




The fluid pressures and free surface displacements are merged: 



I<^y] 


Steps 11 - 13: The dependent displacements for the structure points are 

recovered with Module SDRl, the fluid displacements and pressures 
are merged with the displacement vectors, and the normal MASTRAM 
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output processing is executed. The displacements of fluid free sur- 
face points are printed and plotted. The pressures on the Interior 
fluid points may be printed. All structure displacement, force, 
stress, etc. data Is processed with the standard *"'utlnes. 
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4.2 FUNCTIONAL MODULE FLBMG (FLUID BOUNDARY MATRIX GENERATOR) 


Entry Point ; FLBMG 
Purpose 

To assemble the hydroelastic area factor and gravity stiffness matrices. 
In addition, the hydroelastic USET vector is constructed for defining the 
various structural and fluid sets. 

DMAP Calling Sequence 

FLBMG GE0M2 , ECT , BGPDT , S IL , MPT , GE0M3 , CSTM , USET , EQEXitl/ USETF , US* TS , AF , 

DKGG/S ,N,N0GRAV/S ,N,N0FREE/S,N,TILT $ 

Input Data Blocks 

GE0M2 - Element Connection Data 

ECT - Element Connection Table 

BGPDT - Basic Grid Point Definition Table 

SlL - Scalar Index List 

MPT - Material Properties Table 

GE0M3 - Gravity Load Data 

CSTM - Coordinate System Transformation Matrices 
USET - Displacement Set Definition Table 

EQEXIN - Equivalence between external and internal grid points 
Notes 

1. The CSTM may be purged. 

2. GE0M3 may be purged only if no gravity effects are computed. 

Output Data Blocks 

USETF - Set Definition Table, fluid and structural points 
USETS - Set Definition Table, structural points only 
AF - Fluid Area Factor Matrix 
DKGG - Gravity Stiffness Matrix 


Note 

1. DKGG will be purged if no gravity effect is computed. 
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Parameters 


N0GRAV • Input - integer - no default. Flag which specifies whether 
gravity effects are to be computed (-1 ignores gravity) . 

N0FREE - Output - integer > no default. Flag which specifies whether a 
fluid free surface exists (-1 Implies no free surface) . 

TILT “ Output - complex - default® (1.0, 0.0). Components of the normal 
vector to the free surface. Used as input to SDR2 for plotting. 

Method 

Subroutine FLBMG is a small driver which calls the following routines to 
perform the various phases of matrix generation and assembly. 

1. FLBELM processes the fluid/structural boundary elements and the fluid 
free surface elements. In addition, vectors describing the connec- 
tivity of fluid and structural elements are created for use when 
allocating core during matrix assembly. 

2. FLBSET produces the two new fluid USET vectors for use in the hyc^o- 
elastlc analysis. 

3. FLBEMG calculates the individual element area factor and element 
gravity stiffness matrices. 

4. FLBEMA assembles the element matrices into full size output matrices. 
Subroutines 

The utility subroutines PRETRl) and PREMAT are called by FLBEMG so that 
lower level subroutines may call entry points TR.\NSD and MAT to obtain coordi- 
nate system transformation matrices (CSTM) data and material properties data, 
respectively. The utility routines GMMATD, DCR0SS, and ON0RM are also used 
for basic in-core matrix computations. 

lor purposes of communication between subroutines of the FLBMG module, the 
following common blocks are used. 


Block 


Variables in Order 

/FLBPTR/ 

ERROR 

- Fatal error flag 


IC0RE 

• Next available word of core 


LC0RE 

- Length of open core 


IBGPDT 

- Start of BGPDT data 
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/FLBZZl/ 

/FLBZZ2/ 

/FIEF’LL/ 


NBGPDT - Length of data 
ISIL - Start of SIL data 

NSTL - Length of data 

IGRAV - Start of GRAV data 
NGRAV - Length of data 

IGRID - Start of Grid Point Connectivity Table 

NGRID - Length of table 

IBUFl - GIN0 buffer 1 

IBUF2 - GIN0 buffer 2 

IBUF3 - GIN0 buffer 3 

IBUF4 - GIN0 buffer 4 

IBUF5 - GIN0 buffer 5 


Z - Open core for the first three phases of matrix 

generation. 


Z - Open core for the last phase, or matrix 


assembly. 


GE0M2 

£CT 

BGPDT 

SIL 

MPT 

GE0M3 

CSTM 

USET 

EZEXIN 

USETF 

USETS 

AF 

DKGG 

FBELM 

FRELM 

C0NECT 

AFMAT 

AFDICT 

KGMAT 

KGDICT^ 


GIN0 file numbers for the '^arious input, output 
and scratch files. 
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Subroutine Name: FT.BBLM 


Entry Point : FLBELM 

Purpose: To process the fluid elements constituting the fluid/structure 

boundary and to determine connectivity between these elements. 

Calling Sequence: CALL FLBELM 

Method: The BGPDT data describing grid points is opened and read into 

open core. GE<^M2 is then opened and the record containing CFLSTR bulk 
data cards Is located. The data is then read and stored below the BGPDT 
data in the Element Boundary Table as follows: 


Word No. 


Description 


1 Fluid element ID 

2 Structural element ID 

3-6 (zero) 

7 Gravity ID 

✓ 

The ECT data block is then read and each 2-D element type is processed. 

As each element is processed, if the structural element is connected to 
any fluid element, the grid points of the structural element are Inserted 
Into the four zero words in the element boundary table record. (If the 
structural element is triangular, only three grids are found.) The CFREE 
cards are then read from GE0M2 and stored in a similar seven-word record 
In core. 


Word No. 

1 

2 

3 

4-6 

7 


Description 

Fluid element ID 
-1 

Element FACE ID 
zero 

Gravity ID 


The ECT is then read again and only fluid elements are processed. For each 
element the following steps are performed. 

If the fluid element is on the free surface (the structural element 
ID will be negative): 
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1. Subroutine FLFACE is called to locate the 3 or 4 fluid grid points 
on the desired face which forms the free surface. 

2. A 7-word record is written on file FRELM which contains the fluid 
element ID, material ID, gravity ID, and the fluid grid points 
defining the desired face. 

If the fluid element is connected to a structural element: 

1. Subroutine FLFACE is called to locate the 3 or 4 fluid grid points 
on the face which forms the boundary. 

2. A 12-word record is written on file FBELM which is similar to 
FRELM except the structural element ID and grid points are also 
included. 

As the above steps are being performed a table of grid point connectivity 

is also created. This table provides three pieces of information: 

1. The maximum number of structural grid points connected to each 
fluid point on the boundary. 

2. The maximum number of structural grid points connected to each 
structural point on the boundary. 

3. The maximum number of fluid points connected to each fluid point 
on the free surface. 

The open core layout for FLBELM is as follows: 


IBGPDT 

lELil 


IVEC 

IBUF3 

IBUF2 

IBUFl 


Grid Point Data 


Fluid element data for boundary 
and free surface 


Grid point connectivity table 
GIN0 buffer for FRELM 
GIN0 buffer for ECT 
GIN0 buffer for GE0M2 and FBELM 


4 words/entry 
7 words/entry 
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Subprogram Name: FLFACE 


Entry Point; FLFACE 

Purpose: To locate the fluid grid points which describe the face of a 

fluid element. The particular face may be specified by either a face 
number or the grids of a structural element which form a boundary with 
the face. 

Calling Sequence: CALL FLFACE (TYPE,ECT,EBT,ELT) 

TYPE - The element type 

ECT “ The ECT record for the desired fluid element 

ELT - The element boundary table entry for the fluid element. 

The 7-word block, as described in the previous section, 
contains either the fluid face number or the structural 
grid ID^s. 

GRID - Fluid grid points located (3 or 4) 

Method: If the fact ID is presented: 

The corresponding fluid grid points are obtained directly. The face 
numbering conventions are described in the User’s Manual under the 
bulk data description for each fluid element type. 

If the structural grid points are provided: 

1. Find normal and centroid of structural element face 



R 

c 

-X. 

R 

c 

where 

BGPDT. 


= R2-R1 



for triangle element 
= ^ ^2 ^ ^3 ^ quad element 

are coordinates for structural grid points obtained from the 


Ri2 X 

|Rj ^2 ^ ^13! 


iCRi ^ ''2 
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2. For each face of the fluid element, find the angle between the fluid 
face and structure face. If this angle is less than 30®, also find 
the distance from the centroid of the structure element normal to the 
face of the fluid element. 



^12 

“ ^2 - 

^1 






^^13 



( 


^12 

^^13 



l ’^12 

*'131 


IF Ik, • K I < 0.866 THEN 
' I s' 

Ah = |k^ • (R^ - 

where r^ are coordinates for fluid grid points obtained from the BGPDT. 

3. The face choosen has the smallest Ah value with the angle between 
faces less than 30®. 

Subprogram Name: FLBSET 

Entry Point : FLBSET 

Purpc-^; To construct the hydroelastic USET vectors. 

Calling Sequence: CALL FLBSET 

Method: The SIL list is read into core beneath the BGPDT and the USET 

vector after that. The list of grid point connectivity is then written 
to the CONECT file in two files, one file for each matrix to be assembled 
later. The files consist of three-word records as follows: 

Word No. Description 

1 SIL number 

2 Maximum grid points connected to this SIL 

3 Maximum SIL's connected 

The USET table is then read into core. A list of free surface grid points 
i ; then constructed by examining each entry in the grid point connectivity 
table. The BGPDT is then scanned and the following USET bits are set: 


i - 

»' 
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If fluid point: 

Set Y bit 

Search list of free surface grids; if present, set fr, Z, and a bit, 
otherwise point is interior fluid point and set i bit. 

If structural point: 

Set X and Z bit 

If a bit is set, set a 

Hydroelastic USET bit positions are as follows: 


Set 


Bit Position 


Description 


The 

the 


U 

a 

U 

X 

U 



25 

9 

8 

7 

6 


a + fr 

Structure only 
Fluid only 
Free surface 
X + fr 


z 

U- 5a bits (structure only) 

a 

4 Interior fluid point 

new USET vector is then written to the USETf file. In addition, only 
structure points, x bit set, are written to the USETS file. 


The open core layout for FLBSET is as follows : 


IBGPDT 

ISIL 

lUSET 

IFREE 

IVEC 

IBUFl 


Grid point data 


SIL list 


USET vector 


List of free surface grids 


Grid point connective table 


GIN0 buffer 


4 words/entry 


4-15 




Subprogram Name: FLBPRT 


Entry Point; FLBPRT 
Purpose; 

1. Prints, at user request via DIAG 23, a list of degrees of freedom. 

For each degree of freedom, an indication is made identifying the sets 
to which it belongs. 

2. Prints, at user request via DIAG 24, the contents of selected dis- 
placement sets. For each set, a list of all degrees of freedom 
belonging to the set given. 

Calling Sequence: CALL FLBPRT (IUSET,IC0RE,IBUF) 

lUSKT - Open core location of USETF table 
ICORE - Start of open core for storage of EQEXIN 
IBUF - G1N0 buffer location 

Method: The DIAG flags are tested and local variables set. Table EQEXIN 
is then read into open core and sorted. If DIAG 23 is set. Table USETF 
(already in open core) is examined and the external degree of freedom is 
extracted from EQEXIN and printed along with the set indications. If 
DIAG 24 is set, the transpose process takes place. 

Subprogram Name : FLBEMG 

Entry Point: FLBEMG 

Purpose: To calculate element gravity stiffness and area factor matrices. 

These element matrices, along with dictionary entries to describe them, 
are written to flies for use in the matrix assembly phase. 

Calling Sequence; CALL FLBEMG 

Method: The material property data is read into core below the SIL list 

using utility PREMAT. PRETRD is then called to set up coordinate trans- 
formation data in core. GE0M3 is then opened and GRAV bulk data cards 
are located and read into open cote. The open core layout for FLBEMG will 
then be: 
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IBGPDT 


4 words /entry 


ISIL 

IMPT 

ICSTM 

IGRAV 


IBUF5 

IBUF4 

IBUF3 

IBUF2 

IBUFl 


Grid point data 


SIL list 


Material data 


Transformation data 


Gravity data 


GIN0 buffer for KGDICT 
GIN0 buffer for AFDICT 
GIN0 buffer for KGMAT 
GIN0 buffer for AFMAT 
GIN0 buffer for FBELM 


File FBELM is then read and, for each element, subroutine B0UND is 
called to calculate the area factor matrix and gravity stiffness matrix 
for the element. The gravity stiffness calculations will be performed only 
if gravity loads are present. These matrices are then written to files 
AFMAT and KGMAT respectively. An AFMAT entry contains several 3x1 matrix 
partitions as follows: 


Word No. 

1-3 

4-6 

for triangle faces 
7-60 

for quad faces 
7-102 


Description 

fluid SILS 
structure SILS 

9x3 matrix 

12 X 4 matrix 


The KGMAT file is written in a similar manner except is consists of 
3x3 partitions. Dictionary entries are then written to the AFDICT and 
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KGDICT files to describe each column of the above matrix partitions. 

The file position is also recorded so the matrix partitions may be accessed 
in a random sequence during assembly. A AFDICT entry is as follows: 

Word No. Description 

1 Column SIL number 

2 File position 

After the FBELM file is exhausted, a similar procedure is performed 
with the FRELM fl].e. The major difference is that subroutine FLFREE is 
called to compute the free surface effects on these elements. 


Subprogram Name : FLFREE 

Entry Point: FLFREE 

Purpose: To compute the area factor and gravitational stiffness matrices 

for a fluid element on the free surface. 

Calling Sequence: CALL FLFREE (FRREC,AFE,NAFE,KGE,NKGE) 

FRREC - FRELM record for a fluid element 

AFE - Element area factor matrix. Maximum size is 4 x 4 double 
precision words 

NAFE - Number of words in the AFE matrix 

KGE - Element gravity stiffness matrix. Maximum size is 4 x 4 double 
precision words 

NKGE - Number of words in the KGE matrix 
Method : 

1. If the fluid face is quadratic, it is divided into four overlapping 
triangles for steps 2 through 4. 

2. The area. A, of the fluid triangle is computed. 

3. The area factor terms for this triangle are computed and inserted 
into the full size element matrix. 


lAijl 


for 


A/12 i i j 
A/ 6 i * j 


i - 1,2,3 
j “ 1,2,3 
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4. If gravity loads are requested, the additional stiffness terms are 
computed and inserted Into the element matrix. 

(K^jl - Pg[A^j] 

for 

1 = 1,2,3 
j - 1,2,3 

where 

p * fluid density 
g =» gravitation constant 

5. If the fluid face Is quadratic, the area factor and gravity matrices 
are divided by two to account for the overlapping triangles. 

Subprogram Name : FLBEMA 

Entry Point : FLBEMA 

Purpose: To assemble the element matrices into either the AF or DKGG 

matrix. 

Calling Sequence: CALL FLBEMA (TYPE) 

TYPE - 1 for AF matrix 

2 for DKGG matrix 

Method: One G-size vector is allocated at the top of open core. This 

vector will be the column list, where each entry points to the open core 
location where data for that column is held. The column vector is then 
zeroed and the C0NECT file is opened and positioned for the type of matrix 
being built. As each entry from the C0NECT file is read, the necessary 
core is allocated to hold the terms for that column. A pointer to this 
core is inserted into the column pointer vector. The amount of core 
required for each column is: 

1 + MAX GRIDS + 2 * MAX SILS 

This area will then hold one SIL value and up to 3 double precision matrix 
terms for each grid point entry in the column. 
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This process is continued until the C0HECT file is exhausted or insuf- 
fient core remains to held the next column. The appropriate dictionary 
file is then read and for each entry it is determined if the desired column 
is In core for this pass. If it is, the matrix terms are obtained from the 
matrix term file, xMAT, and stored In core. When all dictionary entries are 
processed, each column in core is written to the output matrix with subroutine 
PAKC0L. Note that only the columns with nonzero terms are stored In core, 
making It necessary to pack out null columns for the others. 

If any non-zero columns remain, l.e., entries remain in the C0NECT file, 
core is reallocated and additional passes are made through dictionary and 
matrix term files to build these columns. 

The open core layout for the assembly phase Is as follows: 


IC0L 

Column pointer 








- Next available space 






+ SIL 


Storage for matrix columns 




3 matrix terms 






- SIL 






1 matrix term 




- 

■V 

BUF2 

GIN0 buffer for AF , DKGG and xMAT 



BUFl 

GIN0 buffer for C0NECT and xDICT 




Subroutine Name: PAKC0L 

Entry Point : PAKC0L 

Purpose: To build the single matrix columns and pack it out to the 

output matrix. 

Calling Sequence: CALL PAKC0L (REC,LREC) 

REC - In-core record containing the SIL's and matrix terms for the 
column 

LREC - Length of this record 


4-20 



Method : 


1. The entries for each SIL (one SIL and up to three matrix terms) are 
sorted on the absolute value of the SIL. 

2. If any duplicate STL's exist, their terms are added together. 

3. The individual matrix terms are then packed to the file using ZBLPKI. 

Subroutine Name: B0UND 

Entry Point: B0UND 

Purpose: To obtain area factor and stiffness matrices as function of 

area common between specified fluid and structural element faces. 

Common Sequence : CALL B0UND (FBREC ,AFE ,NAFE ,KGE,NKGE) 

FBREC(12) - Fluid boundary record: 

(1) Fluid element ID 

(2) Structural element ID 

(3-6) (internal) Grid nos. for structural element face 

(7) Gravity ID 

(8) Material ID 

(9-12) (internal) Grid nos. for fluid element face 
AFE(48) - Area factor metric storage 

NAPE - Number of (single precision) storage locations in AFE 

returned by subroutine B0UND. 

KGE - Stiffness matrix storage. 

NKGE - Number of (single precision) storage locations in KGE 

returned by subroutine B0UND. 

Method: After transforming coordinates of both element faces into 

system in which x/y coordinates of point 1 and y coordinates of point 2 
of the fluid element are zero, the area common to both faces when projected 
on the x,y plane is obtained by determining the coordinates bounding this 
polygon of common area* The load distribution factors and the area 
factors are then calculated by isoparametric interpolation. After any 
required gravitational stiffness terms have been added to the stiffness 
matrix, both the area factor (AFE) and stiffness (KGE) matrices are 
transformed to global coordinate system (if required) and the stiffness 
matrix is rearranged to column storage to conform to NASTRAN requirements. 

Design Requireraents/Limitations : Neither element (face) may have more 

than four edges. If there is no common area, User Warning Message 8014 is 
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printed. If an undefined gravity ID Is referenced, User Fatal Message 8013 
Is printed. I'Hien element geometry Is such that isoparametric interpolation 
falls. User Fatal Message 8005 is printed. 

Subroutine PTINTR (Point of INTe Rsectlon) 

Entry Point : hone 

Purpose: Returns (double precision) xy coordinates of point of inter- 

section (if any) of two lines. 

Calling Sequence; CALL PTINTR(A,AA,B,BB,S,K,EPS) 

A(2) ,AA(2) - xy coordinates of end points of line segments A and AA, 
B(2),BB(2) and B and BB (real, double precision) 

S(2) - xy coordinates (returned) of point of intersection of 

line segments AAA and BBB; S ■ (0,0) if no point of 
intersection (real, double precision). 

K - Condition flag (integer) returned (single precision): 

1, point of intersection is S 
■■0, point of intersection is S, an end point of 
one of the lines 
■ -1, no point of intersection 

EPS (2) - Calculation significance tolerance (real, double 

precision): EPS(l) for area, angle calculations; EPS(2) 

for length calculations. 

X = 0 if |X| 5 € 

Method: K > 0 returned only when point of intersection, S, exists; e.g., 
non-parallel line segments do intersect at S or segments of same parallel 
line have some segment in common (or are within £2 of each other). K • 0 
is returned when distance between S and any end point of AAA or BBB is < C 2 . 

Design Requirements* All lines are defined in two dimensions only; i.e., 
only (x,y) values are considered; all real variables supplied to and returned 
from PTINTR, as well as all internal calculations, are double precision. 

Subroutine L0CPT ( LOCa te PoinT) 

Entry Point : None 

Purpose: Determines location of each of N points, P, relative to surface 

bounded by M points, S. 
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Calling Sequence: CALL L0CPT(N,P,M,S,K,KS,EPS,L0C) 

N “ Number (< of points to locate (integer) 

~ Coordinates of points to locate (real, double precision) 
m - Number of points (< 4) describing boundary of surface. 

S(3,4) - Coordinates of surface bounding points (real, double 
precision) 

K(2,4) - Indices (within S) of end point of surface odges 
KS(3) - Unit vector normal to surface (real, double precision) 

EPS(2) - Calculation significance tolerance (real, double precision): 
EPS(l) for area, angle calculation; EPS (2) for length 
calculations. 

X = 0 if |xl < e 

L0C(4) - Location Indication for each point tested (Integer): 

■ 1, point inside surface boundary 
“ 0, point on surface boundary 
■ -1, point outside surface boundary 

Method: VP, the vector from leading end point o£ each surface edge to 

point to locate is tested against the vector along surface edge. 

The point is within surface boundary (L0C ~ 1) when it is not 'outside* 
any edge or 'on* any edge. 

The point is on a surface edge (L0C ~ 0) when: 

|VP| < C2 

The point is outside a surface edge (L0C -1) when: 

(VE X VP) • k < -Ej or 
(VE • VP) < ^2 
jVE| + Ej < iVP( 

Design Requirements/Limlts: All real arguments and computations are 

^ , double precision. Not more than four points may be positioned relative 

' to a surface with not more than four edges (boundeu by not more than 

four points) . 
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Subroutine P0LYPT 


( P0LYg on PoinTs) 


Entry Points: None 


Purpose: Determine coordinates of points describing the polygon of area 

common to a triangle (structural element) and another figure (fluid element) 
of three or four sides. 


Calling Sequence: CALL P0LYFT (L0CT0F,STEDGE,TR,NGRIDF, FLEDGE, FL.L0CF0S, 

EPS,NP0LY,P) 


L0CT0F(3) 

STEDGE(2,3) 

TR(3,3) 

NGRIDF 
FLEDGE (2, 4) 

FL(3,4) 

L0CF0S(4) 

EPS (2) 


NP0LY 

P(2,7) 


Indication of position of triangle point relative to 
other figure (see Subroutine L0CPT) (integer) . 

Indices of end points of triangle edges (in array TR) 
(Integer) . 

Rectangular coordinates of structural element (real, 
double precision) . 

Number of points for fluid element (integer) . 

Indices (in array FL) of endpoints of fluid element 
edges (integer) . 

Rectangular coordinates of fluid element points (real, 
double precision) . 

Indication of position of fluid points relative to 
structural element (Integer) . 

Calculation significance tolerance (real, double pre- 
cision): EPS(l) for area, angle calculation; EPS(2) for 

length calculations. 

X = 0 if |xl < e 

Number of points (< 7) found to describe polygon of 
common area (integer) . 

x,y coordinates of points describing polygon of 
common ^rea (real, double precision). 


ethod: Starting with first triangle edge, successive points of inter- 

section of (structural) triangle edges and fluid element edges are 
obtained (by Subroutine PTTNTR) . 


Design Requirements/Limitations: All real arguments and computations 

are double precision. 
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Function DAP0LY 


(Rouble Precision Area of a P0LYg on) 


Entry Points : None 


Purpose: Obtain area of a polygon. 

Calling Sequence: A = DAP0LY(N,P) 

N < 10 - Number of points describing polygon (integer, single 
precision) . 

P(2,N) - x,y coordinates of points describing polygon (real, 
double precision). 

Method : 

A = - i y dx 


Contribution from side whose endpoints are P^, is: 

= 1/2 [(y^^ + yj) * (x^ - x^)] 

Design Considerations: 

1. Polygon must not have more than 10 sides. 

2. All real arguments and internal computations are double precision. 


Function DVMAG (Double Precision V^ector MAGn itude) 

Entry Points : None 

Purpose: Obtains magnitude of a vector. 

Calling Sequence: X = DVMAG (V, EPS) 


V(3) - i,j,k coefficients of vector, VI (real, double precision) 

EPS - Calculation significance tolerance (real , double precision) : 
EPS(l) for area, angle calculation; EPS (2) for length 
calculations. 


Method : 


X = 0 if ix| < e 


A = V • V 

DVMAG » /T if A > e > 0 
* 0 if A S e 
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Design Requirements; All real arguments and internal computations are 
double precision. 


Subroutine DCR0SS (Double precision CR0S S product) 

Entry Point : DN0RM 

Purpose: Obtain cross-product of two vectors. 

Calling Sequence: CALL DCR0SS(X,Y,Z) 

X(3) ,Y(3) ,Z(3) - i,j,k coefficients of vectors X,Y,Z. 

Method: Z = X x Y 

Design Considerations: All arguments and internal computations are 

double precision. 


Subroutine DN0RM (Double precision N0RMalization of a vector quantity) 


Entry Point : DCR0SS 

Purpose: Normalize a vector to its magnitude. 

Calling Sequence: CALL DN0RM(X,MAG) 

X(3) - i,j,k coefficients of vector X (real, double precision). 
MAG - Magnitude of vector X (real, double precision). 


Method : 


A = 



+ 


4 


MAG 

X 


A, A > 0 


MAG MAG 


i 

MAG 


Design Considerations: 

1. DN0RM returns un-normalized vector quantity, X, when A < 0. 

2. All arguments and internal computations are double precision. 


Diagnostic Message 

Messages 8000 through 8014 may be Issued by this module. 
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4.3 FUNCTIONAL MODULE GFSMA (GENERAL FLUID/ STRUCTURE MATRIX ASSEMBLER) 

Entry Point : GFSMA 

Purpose 

To assemble the stiffness, mass and pressure transformation matrices in 
a hydroelastic analysis. 

DMAP Calling Sequence 

GFSMA AXY , AFRY ,KYY , DKAA , DKRFFR ,KAA ,MAA ,GM , G0 , USET ,USETF , PHIA ,PHIX , 

LAMA/KMAT ,MMAT , GI A ,P0UT ,HC/ V , N ,N0GRAV/V ,N , N0FREE/ V , Y ,KC0MP/ 
V,Y,C0MPTYP/V,N,F0RM/V,Y,LM0DES $ 

Input Data Blocks 

AXY - Structure/fluid area matrix 

AFRY - Free surface area matrix 

KYY - Fluid stiffness matrix 

DKAA - Structure gravity stiffness matrix 

DKFRFR - Free surface gravity stiffness matrix 

KAA - Reduced structure stiffness matrix >. 

MAA - Reduced structure mass matrix I 

GM - Multipoint constraint transformation matrix > Direct only 

G0 - Omit point transformation matrix 1 

USET - Structure only set definition table ' 

USETF - Fluid and structure set definition table 

PHIA - Solution eigenvectors, A-set 

PHIX - Solution eigenvectors, X-set V Modal only 

LAMA - Solution eigenvalue table ' 

Notes 

1. AFRY and DKFRFR may be purged if no free surface points exist. 

2. DKFRFR and DKAA may be purged if no gravity exists. 

3. GM may be purged if no multipoint constraints exist. 

4. G0 may be purged if no omit points exist. 

5. PHIA, PHIX, and LAMA may be purged if the direct formulation is used. 

6. KAA, MAA, GM. G0, and USET may be purged if the modal formulation is used. 
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Output Data Blocks 


KMAT - Combinaciun fluid/structure stiffness matrix 

MMAT - Combination fluid/structure mass matrix 

GIA - Pressure transformation matrix 

P0UT - Partitioning vector for the modal displacements 

HC - Constraint transformation matrix for incompressible calculations. 
Notes 

1. P0UT will be purged if the direct formulation is used. 

2. HC will be purged if the inconpr^ * sible calculations are not used. 

Parameters 

N0GRAV - Input - integer - no default. Flag which specifies whether 
gravity effects are present (-1 no gravity). 

N0FREE - Input - integer - no default. Flag which specifies whether 
free surface exists (-1 implies no free surface) . 

KC0M? - Input - real - default = 1.0. Compressibility factor. 

C0MPTYP - Input - integer - default = -1. Type of compressibility 
calculation to be used: 

“1 = structure and free surface are coupled with a spring 
to resist volume change 

+1 == incompressible - constraint equation is generated to 
restrict volume change 

F0RM - Input - integer - default = 01. Type of formulation to be used 
-1 = direct formulation 
+1 = modal formulation 

LM0DES - Input - integer - default = -1. Number of structure modes to 
be used in the modal formulation (-1 indicates all available 
modes are used) . 


Method 

Subroutine GFSMA is a small driver which calls the following routines 
based on the input parameters: 
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1. GFSDIR - assembles the combined stiffness and mass matrices using 
the direct formulation. 

2. GFSM0D - assembles the combined stiffness and mass matrices using 
the modal formulation. 


Subroutines 


Subroutine Name: GFSDIR 

Entry Point : GFSDIR 

Purpose: To assemble the fluid matrices using the direct formulation 

method. 

Calling Sequence: CALL GFSDIR 

Method : 

1. The structure/fluid area matrix is reduced to a set using the following 
steps. 

If MFC’s are present, the MPC points are partitioned out. 


[A 1 => 

xy’ 



[A ] = 

ny 


[A ] + 
ny’ 


m my 


If SPC’s are present, the SPG points are partitioned out. 


[A 1 => 

ny 


Al 

A 

sy 


If omits are present: 


(Afyl 


^ 5 !; 

^oy 


[A 1 = (A ] + [G l'*’[A 1 

ay ’ ay o oy 

Subroutines GFSPTN, SSG2B, and CALCV are utilized to perform these 
calculations. 
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2. If free surface points exist, they are merged with the reduced area 
matrix using subroutine GFSMRG. 


A 


ay 


'fry 


“> [A 1 
wy 


3. If SPC points exist on the fluid they are partitioned out from the area 
and fluid stiffness matrices 

[A 1 => [A 1 A 1 

wy wj ' ws 

and this matrix is transposed using subroutine GFSTRN. 


[A. 1 = [A .] 

JW WJ 


[K 1 => 

yy 


K. . 

31 
K . 
S3 


K 


■js 


ss 


4. If no SPC points exist on the fluid, the first fluid point is reduced out 
to remove potential singularities. 


[K 1 => 

yy 


■=11 1 ■^11 
[*^11 i ''ijJ 


generate [H] by Ny transformation) using subroutine GFSH. 


tH] = i 


-1 (n-1) -1 -1 

-1 -1 (n-1) -1 


[A. ] = [H][A ] 

JW ‘ wy 


For C0MPTYP < 0, the compressibility factor is then generated using 
subroutine GFSC0M. This matrix contains the spring used to restrict 
volume change. 
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(A,> - lA^Kl) 

(K I - KC0MP {A,} lA J 
c c c 


For C0MPTYP > 0 a constraint eauation Is generated to restrict 
volume change. Subroutine GFSHC performs these calculations. 


{A^} = [HHi} 


The largest row in Ac,m, is then chosen. The row of matrix 

is null and all other columns, i, have a 1.0 on the diagonal and 

-A ./a in row m. 
ci cm 


IHJ 


0 i -A^./A^„^ 

I cl cin 


5. The pressure transformation matrix is then found by solving the 
following equation: 


■ ‘V 


Matrix utilities FACT0R and SSG3A are used to perform this solution. 

6. If gravity exists, the additional stiffness is added in 

f*aa' " '■'aa'+l^aal 

otherwise, 

[K^J = tK 1 

aa aa 


7. The free surface stiffness is then merged in with subroutine GFSMRG 
if it exists. 



8. The combined stiffness and mass matrices are computed. 


If no SPC’s existed on the fluid for C0MPTYP < 0 
[K ] = tic ] + [K ] 

'■ ww ww c 

for C0MPTYP > 0 

WW c ww c 

[MT] = 

a 1.0 is then added to row m column m of [MT] using subroutine GFSMT. 


ww 


= [MTl + 


1.0 

~o“ 


. 0 

.... 


otherwise, 

[K ] = [K ], [M ] » [MT] 

ww ww ’ ww '■ 

and [K ] is written on data block KMAT and [M ] is written on data 
WW ww 

block MMAT. 

9. The final pressure transformation matrix is created by expanding the j 
size matrix to y size and then removing free surface points. 

If no SPC's exist on the fluid, 

tV’ - 

otherwise merge zeros in. 




1 
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finally partition out the free surface 


(G 1 “O 

‘ yw 


frw 

‘"‘iw 


[G. ] Is then written on data block GIA. 

IW 


Subroutine Name: GFSM0D 


Entry Point : GFSM0D 

Purpose: To assemble the fluid matrices using the modal formulation 

method . 

Calling Sequence: CALL GFSM0D 

Method: 

1. A dummy uset vector, USETD, is created which contains the following 
bit positions: 

Uj, - modal point, UA + UNZ 

- desired modal point (based on LM0DES parameter) 

Un? - modal point to be skipped 
Ufj. - free surface point 
Uh -Ufr+U^ 

Files PHIX and USETF and parameter LM0DES are used to set these bits. 

2. The desired modes are partitioned from the PHIX matrix using the 
USETD vector. 


IPHIX) 


10 r ' 0 cl 

xC • xn? 


3. The area matrix is then transf^rmed. 


‘V’ 


4. If the free surface points exist, they are merged in using subroutine 
GFSMRG. 
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Ar 

-iX. 

^fry 




5. If SPC points exist on the fluid they are partitioned out from the 
area and fluid stiffness matrix. 


‘V 


and this matrix is transposed using subroutine GFSTRN. 

,T 




‘V) 


[K ] 

yy 


ii I js 


K 


si ss 


6. If no SPC points exist on the fluid, the first fluid point is reduced 
out to remove potential singularities. 


[K 



K 

K 


11 

il 


K 

K 


li 

ii 


The [H] matrix is then generated using subroutine GFSM. 
Ny transformation) 


(N 


(y-1) 


by 


[H] 


-1 (n-1) 



-1 -1 . . . 

(n™!) **jL . . . 




The compressibility matrix which contains the spring is generated to 
restrict volume change. 


{Ac> = [A^yHl} 

[K^] = KC0MP {A^}lA^J 
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7. The pressure transform matrix is then found by solving the following 
equation: 

IKjjltOj,,] - (AjJ 

Matrix utilities FACT0R and SSG3A are used to perform this 
solution. 

8. The generalized stiffness and mass terms are tb' extracted from 
the LAMA data block to form the diagonal modal mass and stiffness 
matrices. 


|K„1 


h 






M 


1 . 


M 


9. If free surface points exist, the modal mass matrix Is expanded and 
the final mass matrix is formed. 


1° 

0 ! 0 
I j 


=> “ modes + free surface) 


‘V - i\h' * ‘V 


and is written on data block MMAT. 

10. If gravtly exists, the additional stiffness Is transformed and added In. 

tPHIA] => «*an^J 


[K^J « [K_J + (0„rl IDK J[0 


'a 


a 
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otherwise 


“'u’ ■ 


11. The free surface stiffness is then merged in if it exists. 


K-- I 



I DK, . 

I frfr 

I 




12. The final stiffness matrix is then ccmpuced by adding in the 
compressibility term if it exists. 

If no SPC's existed on the fluid 


otherwise 




and 


i\hi 


is written on data block KMAT. 


13. The final pressure transformation matrix is created by expanding 
the J size to h and then removing the free surface points. 

If no SPC's exist on the fluid. 




otherwise zeros are merged in. 


fit 


■> «=yhl 


Finally the free surface is removed. 




frh 

^Th' 


[C^j^] is written on data block GIA. 
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Subroutine Name: GFSSPC 


Entry Point: CFSSPC (in GFSUTL) 

Purpose: To calculate a partitioning vector to remove the first row 

and columns of the fluid stiffness matrix if no SPC's are provided. 

Calling Sequence: CALL GFSSPC (NUT, PVEC) 

NUY ~ Size of fluid stiffness matrix 

PVEC - Gino file number of partitioning vector 

Subroutine Name : GFSC0M 

Entry Point: GFSC0M (in GFSUTL) 

Purpose: To compute the fluii comprassibility matrix fo* the .ase 

where no SPC's exist on the fluid. This matrix contains the sprang 
factors that couple the free surface and the structure points to resist 
volume changes. 

Calling Sequence: CALL GFSC0M(AWY,NUY,KC,1DENT,AC,SCR) 

AWY - Gino file number of whe ANY matrix, direct formulation, or 
AMY matrix, modal formulation 
NUY - Size of fluid matrices 

KC - Gino fixe number of the compressibility matrix to be 
calculated 

IDENT - G^no file number of a scratch file 

AC - Gino file number of a scratch file 

SCR - Gino file n’xmber of a scratch file 

Method: See Step 4 under subroutine GFSDIR for the calculations used 

to geneiate KC. 

Subroutine Name: GFSH 

Entry Point: GFSH (in GFSUTL) 

Purpose: To calculate the H transformation matrix used when the first 

row and column of the fluid stiffness matrix is removed. 

Calling Sequence: CALL GFSH(NlfY,H) 
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NUY - Size of fluid stiffness matrix 
H - Glno file number for H matrix 

Method; See Step 4 under subroutine GFSDIR for the calculations used 
to generate H. 

Subroutine Name : GFSHC 

Entry Point: GFSHC (in GFSUTL) 

Purpose: To generate the constraint transformation matrix for the pure 

incompressible calculations when requested by the user. These constraints 
will restrict any volume changes by the fluid. 

Calling Sequence; CALL GFSHC(AWY,NUY,HC,IDENT,AC,MR0W) 

AWY - Gino file number of the AWY matrix 

NUY - Size of fluid stiffness matrix 

HC - Gino file number of the HC matrix 

IDENT - Gino file number of a scratch matrix 

AC - Gino file number of a scratch matrix 

MR0W - Row number of the largest term in the calculated AC matrix 

Method: See Step 4 under subroutine GFSDIR for the calculations used to 

generate HC. 

Subroutine Name; GFSMT 

Entry Point: GFSMT (in GFSUTL) 

Purpose: To add a 1.0 value to the null row and column of the fluid 

structure mass matrix to prevent singularities. The routine is only used 
when the incompressible calculations are requested. 

Calling Sequence: CALL GFSMT (MT,MMAT,MR0W) 

MT - Gino file number of the input mass matrix 
MMAT - Gino file number of the final mass matrix 
MR0W - Null row and column number. Output from subroutine GFSHC. 

Method: Subroutine CPYSTR is used to copy each column up to column MR0W 

from matrix MT to MMAT. A 1.0 is then packed in matrix MMAT in row 
position MR0W using subroutine ZBLPKI. The remainder of matrix MT is 
then copied to MMAT usine subroutine CPYFIL. 



Subroutine Name: GFSPTN 


Entry Point : GFSPTN 

Purpose: General purpose partition routine to perform the following 

partition: 


(A] -=> 


I 


All A12 


A21 1 A22 


Calling Sequence: CALL GFSPTN (A, All , A21 ,A12,A22, RPART, CPARY) 

A - Gino file number of matrix to be partitioned 

All ) 

A2 X f 

> - Gino file numbers of output matrices. Any of there may be 
A22 ) zero if that partition is not desired. 

RPART - A partitioning vector whose length is the size of a row in 

A. It is used to partition the columns of A and may be zero 

if A12 and A22 do not exist. 

CPART - A partitioning vector whose length is the size of a column in 

A. It is used to partition the rows of A and may be zero if 

A21 and A22 do not exist. 


Method: The input arguments are used to initialize common block /PARMEG/ 

and subroutine PARTN is called to perform the actual partition operation. 


Subroutine Name: GFSMRG 


En try Point: GFSMRG 

Purpose: General purpose merge routine to perform the following merge: 


All ! A12 

I 

A21 I A22 


-> [A] 


Calling Sequence: CA1.L GFSMRG(A, All, A21,A12,A22, RPART, CPART) 

Argument values are the same as those discussed under subroutine 
GFSPTN. 



< 


Method: The input parameters are used to initialize common block /PARMEG/ 

and subroutine MERGE is called t "erfonn the actual merge operation. 

Subroutine Name: GFSTRN 

Entry Point ; GFSTRN 

Purpose; To transpose a matrix. 

Calling Sequence: CALL GFSTRN(A,AT,I,SCR) 

A ~ Gino file number of matrix to be transposed 
AT - Gino file number of the transpose 
I - Gino file number of a scratch file 
SCR - Gino file number of a scratch file 

Method: Subroutine SSG2B is used to find the transpose by solving the 

following equation: 

[AT] = [A]'^[I] 

^diere I is an identity matrix. 

This procedure is more efficient at finding the transpose then subroutine 
TRANPl when matrix A is sparse. 

Subroutine Name: GFSWCH 


Entry Point: GFSWCH 


Purpose: To switch the units on which two files reside. 


Calling Sequence: CALL GFSWCH (FILE1,FILE2) 


FILEl 

FILE2 


Gino file numbers of files to be switched. 


Method: The switch is performed by updating the FIST and FIAT to point 

to the new files. A check is made to ensure that any stacked data blocks 
in the FIAT are also changed to reflect the switch. 


Diagnostic Messages 


A number of messages produced by the matrix utilities can be issued by 
this module. 
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4.4 UTILITY MODULE 'TRAILER' 

Name : TRAILER 

Purpose 

To examine or modify the trailer of a GIN0 data block. 

DMAP Calling Sequence 

TRAILER A//C,N,opt/C,N,word/5,N,value $ 

Input Data Blocks 

A - Data block for which the trailer is desired 


Output Data Blocks : None 


Parameters 

1. opt is a BCD operation code from the list below, (input, no default) 

RETURN - The value of the specified trailer word is to be returned 
ST0RE - The value of the specified trailer word is to be changed 

2. word is an integer value which specifies the particular trailer word 
to be operated on. The value of this parameter must be 1 < word ^ 6. 
(input, no default) 

3. value is an integer parameter which will contain the value of specified 
trailer word, opt = RETURN, or contains the value to be stored in the 
trailer, opt = ST0RE. If specified, the range must be 0 ^ value < 65535. 
(input, no default) 


Remarks 

1. For matrix data blocks, the trailer positions contain: 


WORD 1 - Number of columns 

WORD 2 - Number of rows 

WORD 3 - Matrix form 

WORD 4 - Type of matrix elements 

WORD 5 - Maximum number of non-zero words in any one column 
WORD 6 - Matrix density x 10^ 


$ 
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2. If the data block is purged, the parameter value will be returned 
negative. 


Examples 

TRAILER M1//C,Y,RETURN/C,Y,1/V,Y,C0L $ 
SAVE C0L $ 

TRAILER M2//C,Y,ST0RE/C,Y,3/C,Y,1 $ 
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4.5 DMi\P EXECUTIVE OPERATION MODULE C0MP0N 


Name ; C0ND0N (compilation on) 

Purpose 

To allow blocks of DMAP statements to be compiled or skipped depending 
on bulk data parameter value. 

DMAP Calling Sequence 

C0MP0N n,param $ 

C0MP0N Cjparam $ 

where : 

1. n is a BCD name of a label which specifies the end of the DMAP 
statement block. 

2. c is an integer constant which specifies the number of DMAP statements 
in the block. 

3. param is the name of parameter that appears on a PARAM bulk data card. 


Example 


C0MP0N 

END, PI $ 

M0DULE1 

A/B/V,Y,PN $ 

• 

M0DULEN 

A/B/V,Y,PN $ 

LABEL 

END $ 

C0MP0N 

2, PI $ 

M0DULE1 

A/B/V,Y,PN $ 

M0DULE2 

A/B/V,Y,PN $ 


Remarks 

1. The block of DMAP statements specified by the label or count is 
skipped if the value of parameter is false (param > 0) . If the 
parameter value is true (param < 0) , the block of DMAP statements 
will be compiled. 

2. If no PARAM bulk data card is provided, a value of false is assumed. 
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3. If the form of C0MP0N specifying a label is used, the label may not 
be specified by any other DMAP instructions including other C0MP0N or 
C0MP0PF instructions. 

4. Comment cards are not included in the statement count. 

5. C0MP0N and C0MP0FF instructions may be nested up to five levels using 
the same rules as for a F0RTRAN D0 loop. 
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4.6 DMAP EXECUTIVE OPERATION MODULE C0MP0FF 


Name ; C0MP0FF (compilation off) 

Purpose 

To allow blocks of DMAP st.:tements to be compiled or skipped depending 
on a bulk data parameter value. 

DMAP Calling Sequence 

C0MF0FF n,param $ 

C0MP0FF c,param $ 

where : 

1. n is a BCD name of a label which specifies the end of the DMAP 
statement block. 

2. c is an integer constant which specifies the number of DMAP 
statements in the block. 

3. param is the name of parameter that appears on a PARAM bulk data 
card. 


Example 


C0MP0FF 

END, PI $ 

M0DULE1 

• 

A/B/V,Y,PN $ 

• 

M0DULEN 

A/B/V,Y,PN $ 

LABEL 

END $ 

C0MP0FF 

2, PI $ 

M0DULE1 

A/B/V,Y,PN $ 

M0DULE2 

A/B/V,Y,PN $ 


Remarks 

1. The block of DMA? statements specified by the label or count is 
skipped if the value of the parameter is true (param < 0) . If the 
parameter value is false (param > 0) , the block of DMAP statements 
will be compiled. 
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2. If no FARAM bulk data card is provided, a value of false Is assumed. 

3. If the form of C0MP0FF specifying a label is used, the label may not 
be specified by any other DMAP instructions including other C0MP0N 
or C0MP0FF instructions. 

4. Comment cards are not included in the statement count. 

5. C0MP0N and C0MP0FF Instructions may be nested up to five levels 
using the same rules as for a F0RTRAN D0 loop. 
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4.7 MODIFICATIONS TO EXISTING NASTRAN SUBROUTINES 


Changes to the standard Level 16 were required to perform the following 
operations: 

1. The three-dimensional fluid elements are connected to special fluid 
grid points with one degree of freedom. Implementing these elements 
required changes to the IFF, GPl, TAl, PL0T, and EMG modules. The 
fluid material definition required changes to the PREMAT utility 
subroutine. 

2. Generating plots of the free surface displacements required changes 
to the SDR2 module to convert scalar displacements to vectors normal 
to the free su.face. 

3. The two-step eigenvalue processing In the modal formulation required 
changes to the CASE module to process separate subcases for eigenvalue 
methods. 

4. The 'GIVENS' option for eigenvalue extraction in the READ module was 
modified to provide more efficient processing and to correct for 
existing errors. Specifically, the double precision operations were 
corrected to provide consistent calculations and the 'spill logic' 
in the trldlagonalizaclon procedure was completely recoded to 
eliminate unnecessary core transfers and I/O operations. 
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FUNCTIOtiAL MODULE CASE (SIMPLIFY CASE CONTROL) 


^4.56 FUNCTIONAL MODULE CASE (SlHPuIFY CASE CONTROL) 

4.56.1 Entry Point ; CASE 

4.56.2 Purpose 

To remove looping considerations from later dynamics modules. 

4.56.3 DMAP Calling Sequence 

CASE CASECC,PSOL/CASEXX/C,N,APPR{JACH/V,N,REPEAT/V.N,L00P $ 

4.56.4 Input Data Blocks 

CASECC > Case Control Data Table. 

PSOL - Power Spectral Density List. 

Note ; PSDL is used only if APPROACH = FREQRESP and Random Analysis is selected in CASECC. 

4.56.5 Output Data Blocks 

CASEXX - Case Control data table for dynamics problems. 

Note ; CASEXX cannot be pure ’ 

4.56.6 Parameters 

APPROACH - Input-BCO-no default. Defines the approach to be used for looping 
criteria. 


BCD Value 

LggP 

STATICS 

N0NE 

1 

REIGEN 


OSO 

N0ME 

0S1 

NONE 

FREQRESP 

DIRECT INPUT MATRICES OR TRANSFER FUNCTIONS 

TRANRESP 

L0AOS 

BLKO 

NONE 

3LK1 

NONE 

CEIGEN 

DIRECT INPUT MATRICES OR TRANSFER FUNCTIONS 
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BCD Value LOOP 

PLA N0NE 

REPEAT - Input and output-InteseT'Set equal to zero outside of the DMAP looi, by the PARAH 
module. -1 If no additi-'iel loops; + loop court if loops. 

L00P “ Output-integer-default » -1. -1 if this is not a looping problem, 0 if this 
is a looping problem. 

4.56.7 Method 

The method of operation depends upon the input parameter APPR0ACH. 

4.56.7.1 Transient Response 

If APPR0ACH • TRANRESP, CASECC is skipped over REPEAT records. If REPEAT - 0, REPEAT is set 
to 1. One record of CASECC is read and copied onto CASEXX. An attempt is made to read another 
record. If no more records exist, REPEAT is set to -1. Also, if this is the first entry to CASE 

(i.e., REPEAT » 1), L00P is set to -1. If additional records exist, REPEAT and L00P are set to 1. , 

i,. 

4.56.7.2 Complex Eigenvalue Analysis 

If APPROACH =• CEIGEN, REPEAT records are skipped in CASECC. If REPEAT « 0, REPEAT is set to 
1. One record of CASECC is read and copied onto CASEXX. The names of the Direct Input Matrices 
and Transfer Functions sets are saved. An attempt is made to read another record. If no 
more exist, REPEAT is set to -1. Also if this is the first entry (i.e., REPEAT > 1) L00P is set 
to -1. If additional records exist, their Direct Input Matrices and Transfer Functions sets are 
compared to those saved. If they all agree, this record is copied onto CASEXX and the process 
is repeated. If they do not agree, REPEAT is incremented by 1, L00P is set to 1, and CASE returns. 

4.55.7.3 Frequency Response 

If APPP0ACH » FREQRESP, the method used is tne same as Complex Eigenvalue Analysis except a 
test is also made for frequency set selection changes. In addition, if PANDPS cards are selected, 
the selected set is read from PSOL and the unique subcase "id’s" referenced are stored. Each sub- 
case id copied onto CASEXX is compared to this list, and the entry is marked as found. If at the 
completion of CASE unmarked entries exist, the routine terminates with message 3033. 
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STRUCTURAL ELEMENT DESCRIPTIONS 


The "stress" data recovery for heat transfer analysis Is performed by subroutines ShHTFI, 
SOHTFF, and SDHTF2 of module S0R2. In Phase I, the matrix K and the matrix C are calculated 
where K Is the 3x3 material matrix and Is a 3 by number of points matrix. For the 
tetrahedron: 

”21 “22 ^3 

“31 ”32 ”33 “34 

_»41 «42 «43 »44_ 

For the WEDGE. HEXA1, and HEXA2 elements, the [C] matrix Is calculated for each subelement. 
The column corresponding to each point of the tetrahedron Is added to the column of the matrix 
corresponding to that point In the whole element. The results are divided by the number of 
subelements. 


[Cg3 - 


In Phase 2, the temperature gri'iient vector and flux vector are calculated with the equations: 

{AT} « [C] {u} 

{q> » -[K](AT}, 


where {u} Is the vector of temperatures of the connected points. 


8.17.3 Hydroelastic Calculations for Solid Elements 


The "stiffness" matrix for hydroelastic problems is generated for the FHSXl, 
FHEX2, FTETRA, and FWEDGE elements (which are Identical in format to the solid 
structure equivalents). The matrix is Identical to the heat transfer matrix 
except 



» K 

- K 

‘xx 

yy 

zz 


* K 

- K 

xy 

xz 

y* 


1/P 

0 


The density, p, is obtained from the MAT subroutine with INFLAG* 9. 
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INFLAG ■ 8 — INFLAG ■ 8 Is us^id only by two-dimensional element subroutines 
In (nodules PLA3 and PLA4. The fourth word of /WilIN/, PLAARG (see below). Is 
stress (o) and Is used as the ordinate In an •'nverse interpolation table look- 
up to obtain the abscissa which Is strain (c). 

If either: a) the ordinate Is in the range of the piecewise linear function 
defined by the table on a TA6LES1 bulk data card, or b) the ordinate Is greater 
than the maximum (which Is also the last) ordinate 1n the table but the slope 
of the line segment jolning'the last two points of the table Is nonzero, then 
the second word of /HAT0UT/ Is set to zero and the abscissa, obtained by Inverse 
linear interpolation or extrapolation. Is stored In the first word of /MAT0UT/. 
If either: a) the ordinate Is less than the minimum (which Is also the first) 
ordinate In the table, or b) the ordinate Is greater than the maximum ordinate 
In the table and the slope of the line segment joining the last two points of 
the table Is zero, then the Integer "1" is stored In the second word of /MATS3UT/ 
(and the first word of /MAT0UT/ Is set to zero). Only MATl cards are searched 
to match the Input MATID. 

TEMP - Average element temperature. Used as the independent variable in a table look- 
up when It Is determined that a material property is temperature dependent. 

Not used when INFLAG ■ 3 or 6. 

PLAARG - Element strain. Used as the Independent variable In a table look-up when E, 

the irz)dulus of elasticity. Is defined as the first derivative of a strain-stress 
curve. Used only in the Piecewise Linear Analysis Rigid Format and only by 
modules PLA3 and PLA4. 

SINTH - Sine of the material property orientation angle. Used only when INFLAG ■ 2 

and the MATIO is found among the MAT2 cards. Used to construct the [U] matrix 
referenced above. 

C0STH - Cosine of the material property orientation angle. The comments on SINTH, 
above, also apply here. 

CptW^N/MAT^UT/ - (output Common Block). Length 20 words. Depending upon the values of INFLAG, 
the output common block is defined variously as follows: 

XKFL AS - xs usd t>~:> c\en\eY^e.. 

'TVievaUe \s rdixmeci ‘lO /Hf\T(^QT/. 
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8. Strain Functional Value (INFLAG =» 8) 


Word Symbol 

1 PLAANS 

? ICELL2 

3-26 


Definition 


Value of strain (e) as an inverse function 
of stress (a) 

= 0 if the input stress is in the range of 
the function 

- 1 if the input stress is outside the range 
of the function 

Undefined 


3*4,3C.4 Method 

and Hf^TF . 

1. FREMAT: All the MATl » MAT3^cards are read from the MPT data block into open I 

core so that each card is assigned 1 + 3*N words of core where N» a function of the card type, is 
the number of material property data items on that card type. The first word is the material iden* 
tification number and each material property is allocated 3 words: the fir.st the input material 
property; the second a table (function) number which gives this material property as a function of 
temperature; the third a table number which gives this material property as a function of stress. 
Initially words 2 and 3 are set to zero. Although the third word is currently used only for MAT1 
cards and for E, the modulus of elasticity, on that card, future development may make use of a 
more general application of stress dependent material properties. If there are r temperature 
dependent material properties for a non-Piecewise Linear Analysis problem, PREMAT is wrapped up 
and a RETURN to the calling routine is executed. 

For a non-Piecewise Linear Analysis problem for which a temperature set for material proper- 
ties was selected in the user's Case Control Deck, all MATTl, MATT2 and MATT3 cards are read into 
open core from the MPT data block. For a Piecewise Linear Analysis problem MATSl cards are read 
into open core from the MPT. A sorted list, with duplicates discarded, of the table numbers 
referenced on these cards is constructed in open core. This* table number list is constructed so 
that every referenced table has eleven locations allocated to it. These eleven locations are used 
as a dictionary for the tables. The contents are: the table number (word 1); the table type 1,2,3 
or 4 (v/ord 2); pointers to the first and last entries in the table (words 3 and 4); parameters 
from the TABLE card (words 5 through 11). Ihe DIT data block is then read. For each table read, 
it is determined by scanning the table number list whether or not the table is required for prob- 
lem solution. If it is required, the table is read into oper. core and the dictionary entry for 
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EXECUTIVE PREFACE MODULE IFP 


Table 1(f). Bulk Data Cards Processed by IFP. 


A 

B 

C 

D 

E 

F 

G 

H 

I 

J 

K 

L M 

N 

291 

CFFREE 

8 

GE0M2 

0 

4 

8 

505 

1 

4404 

75 

1 

-1 

292 

CFLSTR 

8 

GEOM2 

0 

-4 

9 

-1 

1 

4504 

76 

3 

-1 

293 

CFHEXl 

8 

GE0M1 

0 

16 

16 

531 

1 

4604 

77 

4 

-1 

294 

CFHEX2 

8 

GEj)M2 

0 

16 

16 

531 

1 

4704 

78 

4 

-1 

295 

CFTETRA 

8 

GE0M2 

0 

8 

8 

337 

1 

4804 

79 

4 

-1 

296 

CR^EDGE 

8 

GE0M2 

0 

8 

8 

525 

1 

4904 

80 

4 

-1 

297 

^LATF 

3 

GE0M2 

0 

4 

8 

198 

1 

5004 

81 

1 

-1 
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FUMCnONAL MODULE GPl (GEOMETRY PROCESSOR - PHASE 1) 


{k} X {V,} 

{j} = IfiJTlvTT ’ vector); (21) 

{<} * {j} X {k} , ("X* unit vector). (22) 


The orientation of the axes is defined by the matrix 




i, S 


’3 ^3 •'3 


(23) 


5. On each pass of the C0ROij data at least one new system must be converted. After 
each pass the referenced GRIO data is checked and converted. The resulting C0ROij data 
will be the CSTM data block with each entry reduced from 16 to 14 words. 


4.21.7.5 Construction of the 8GP0T, the SIL and the Second Logical Record of the EQEXIN. 

The BGPuT and the SIL data blocks are formed simultaneously. The SIL data block is simply 
a list of the first scalar index for each grid or scalar point. The number of scalar indices (or 
degrees of freedom) for each point is determined by examining the elements connected to each point. 

The maximum number of degrees of freedom for each element type is listed in /GPTAl/. The maximum 
degrees of freedom for each point is determined by reading data block GE0M2 and examining the 
connection information in conjunction with the degree of freedom infomation in /GPTAl/. 

The GPDT data are read a point at a time. The basic location coordinates of the point are 
formed using Equation 8 through Equation 14 and these data are written on the BGPDT file. The SIL 
value for the next point is calculated by incr renting the last value by six (grid point) or by 
one (scalar point), CK^VUv<Ag.\evAed’. : 

A test Is inade on the value of the displacement coordinate system (field 6) in the 
If this value is the integer, -1, the point is a special RINGFL, GRIDF, or GRIDS fluid | 
is given one scalar index, the displacement coordinate system is basic (0), and its location 
coordinates in the BGPDT data block are calculated like a normal grid point. 

Finally the second logical record of EQEXIN is written. This record contains pairs of ex- 
ternal numbers, lO^scalar index + type where type * 1 for a grid point, 2 for a scalar point. 
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opened, if not yet opened, and the next vector present is unpacked into core. 

Data Items are now assembled for the identification record, and this identification 
record is output to the output data block. Output line entries for the point-ID*s requested 
are then written on the output data block forning a data record. At this time, if the user 
requested magnitude/phase for complex outputs, the magnitude/phase computations are per- 
formed on the real/imaginary pairs. 

When all requests have been processed for this vector, the next Case Control record 
is read. If no more Case Control records exist and there are more vectors present, those 
vectors are processed using the last Case Control record's specifications. 

When all vectors have been processed for the current loop pass, the next pass may 
be made for forces of single-point constraint or loads. 

If deformed structure plots are requested, an output plot data block is formed during 
the first loop pass, described above, containing translation components of the displacement 
vector rotated to basic coordinates. Sca\df ^ 

“VVe S cure 

.46.8.21 Subroutine Name: SDR20 c^''op\^ce«AO\^. 

1.. Entry Point: SDR2D 

2. Purpose: To perform stage V as defined above, under "Method". 

3. Calling Sequence: CALL SDR2D 

4. Method: SDR2D performs the phase 2 stress and force recovery computations. In this 
phase actual stresses and forces are computed for the user-requested elements. These 
stresses and forces are a function of the stress matrices computed in Stage III and the 
displacements at the grid points of the elements. 

The operations of SDR20 are dependent upon the Rigid Format being executed. In all cases 
the Case Control data block is opened first. For eigenvalue problems a list of eigen- 
values and mode numbers is read into core from LArtA and CLAMA. For Differential Stiffness 
or Buckling phase 1 problems, the first record of Case Control, which is used in phase 0 
of Buckling or Differential Stiffness, is skipped. For frequency or transient response 
problems, a list of frequencies or times is read into core from PPF or PPT, 
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FUNCTIOMAl MODULE READ (REAL EIGENVALUE ANALYSIS - DISPLACEMENT) 

GIN0 file number of the input matrix - integer - input. 

GINO file numbers of scratch files - integer - input. 

'0* indicates no restart is being made - integer - input. 

Order of the problem - integer - output. 

Frequency range for computation of eigenvectors - real - input. 
Eigenvalue sort order flag - integer - input. 

GINSI file name of the associated data blocks - integer - input. 
Number of eigenvectors to compute - integer - input. 

Number of rigid body modes previously found - integer - input. 
Number of fails to converge on eigenvectors - integer - output. 
Number of fails to converge on eigenvalues - integer - output. 
Maximum number of QR iterations allowed - integer - input. 
Reason for termination - integer - input. 

Open core for VALVEC. 

XI, X9, X18, YY, X8 are dummy variables and are not currently used. 

4.48.8.30 Subroutine Name: SMLEIG 

1. Entry Point: SMLEIG 

2. Purpose: To compute the eigenvalues for a 1 by 1 and 2 bv 2 matrix and the eigenvector 
for a 1 by 1 matrix. 

3. Calling Sequence: CALL SMLEIG (0,0, VAL) 

0 - Array of diagonal values - double precision - output, 

jj - Array of off-diagonal values - double precision - output. 

VAL - Array of eigenvalues - double precision - output. 

4.48.8.31 Subroutine TRIOI 
1. Entry Point: TRIOI 

4.48-19a (12-1-59) 


M{3 

MO, MRl, Ml, M2, M3, H4 - 

RSTRT 

N 

LFREQ. HFREQ 
ORDER 

LAfiA, 0EIGS, PHIA 
NV 

NF0UND 

NVER 

NEVER 

MAX 

ITER 

2 
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FUNCTIONAL MODULE READ (REAL EIGENVALUE ANALYSIS - DISPLACEMENT) 


2. Purpose: To tridiagonalize a symmetric matrix. 

3. Calling Sequence: CALL TRIDI (D,0,V,B,A,L0C,DA) 

D - Diagonal terms of the tridiagonal matrix - double precision - output. 

0 - Off-diagonal terms of the tridiagonal matrix - double precision - 

output. 

V - Scratch *;rray which contains another copy of the diagonal terms at 
the conclusion of TRIDI - double precision - output. 

A - Remainder of core - single precision - scratch. 

B - Scratch array which contains the square of the off-diagonal terms - 
double precision - output. 

DA - Remainder of core - double precision - scratch. 

L0C - Remainder of core - integer - scratch. 

4.48.8.32 Subroutine Name: SIC0X 

1. Entry Point: SIC0X 

2. Purpose: To initialize the arrays in SIC0D and SIC0S. See section 

4.48.8.33. 

3. Calling Sequence: CALL SIC0X (SIN,C0S) 

SIN - Array of sine rotation factors - double precision - Input/output. 

C0S - Array of cosine rotation factors - double precision - input/output. 

4.48.8.33 Subroutine Name: SIC0X 

1. Entry Point: SIC0D 

2. Purpose: To compute the rotation factors for a given row in double 

precision. 

3. Calling Sequence: CALL SIC0D (R0W,D,R0T) 

R0W - The number of the current row to rotate - integer - input. 

R0T - If no rotations are requited for this row, R0T ® 0. Otherwise, 

R0T = 1 - Integer - input. 

D - Array of terms for this row - double precision - input/output. 
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FUNCTIONAL MODULE READ (REAL EIGENVALUE ANALYSIS - DISPLACEMENT) 

4.A8.8.34 Subroutine Name: SIC0X 

1. Entry Point: SIC0S 

2. Purpose: To compute the rotation factors for a given row in single 

precision. 

3. Calling Sequence: CALL SIC0S (R0W,DS,R0T) 

R0W - The number of the current row to rotate - Integer - input. 

DS - Array of terms for this row - single precision - input/output. 

R0T - If no rotations are required for this row, R0T=O. Otherwise, 

R0T = 1 - integer. 

4.48.8.35 Subroutine Name: R0TAX 

1. Entry Point: R0TAX 

2. Purpose: To initialize the arrays in R0TATE. See section 4.48.8.35. 

3. Calling Sequence: CALL R0TAX (SIN,C0S) 

SIN - Array of sine rotation factors ~ double precision - input. 

C0S - Array of cosine rotation factors - double precision - input. 

4.48.8.36 Subroutine Name: R0TAX 

1. Entry Point: R0TATD 

2. Purpose: To rotate as much of the matrix as fits into core using 

double precision arithmetic. 

3. Calling Sequence: R0TATD (0,A,R0W,R0W1,R0W2,L0C) 

0 - Array of off-diagonal values - double precision - input/output. 

A - Partition of the matrix held in core - double precision - Input/output. 

R0W - The row number of current rotation row - integer - input. 

R0W1 - The row number of the first row of the matrix partition in core - 
integer - input 

R0W2 - The row number of the last row of the matrix partition in core - 
integer - input. 
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L0C - For spill case a vector containing the open core location for each 
row currently in core - integer - input. 

4.48.8.37 Subroutine Name: R0TAX 

1. Entry Point; R0TATS 

2. Purpose: To rotate as much of the matrix as fits into core using 

single precision arithmetic. 

3. Calling Sequence: CALL R0TATS (0S,AS,R0W,R0W1,R0W2,L0C) 

0S - Array of off-diagonal values - single precision - input/output. 

AS - Partition of the matrix held in core - single precision - input/ 

output. 

R0W - The row number of the current row - integer - input. 

R0W1 - The row number of the first row of the matrix partition in core - 
Integer - input. 

R0W2 - The row number of the last row of the matrix partition In core - 
integer - output. 

L0C - For spill case a vector containing the open core locations for 
each row currently in core - Integer - input. 

4.48.8.38 Subroutine QRITER 

1. Entry Point: QRITER 

2. Purpose: To obtain the eigenvalues of a tridiagonal matrix by the 

Ortega - Kaiser QR iteration technique. 

3. Calling Sequence: CALL QRITER (VAL,0,L0C,QR) 
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INTRODUCTION 


This document Is the Programmer's Manual for the MESHGEN mesh generation 
program. It Includes descriptions of all routines used by MESHGEN as well 
as the data base and file operations. Additionally, a section on system 
considerations points out areas of code that would require change for 
Implementation on other computer hardware. 


1. MESHGEN ORGANIZATION 


MESHGEN has been developed and organized In a modular manner. Such a design 
allows both ease in debugging and simplicity of modification. The main 
driver routine performs virtually no computation, but serves only as a 
logical framework to guide execution. The overall logic of this routine 
is shown In Figure 1 with annotations describing functions. 

For the most part, all data Is passed between major routines through common 
blocks and, in some cases, by scratch files. This again enhances the 
ability to expand and modify the system. 
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FIGURE 1. OVERALL FLOW OF MESHGEN 
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FIGURE 1. OVERALL FLOW OF MESHGEN, Cont’d 
Description of Functions 


PREPAS and DATAEX are used to prepare the input MESHLAN 
sequence for processing. 

LEXER performs the lexical analysis of the MESHLAN sequence 
and builds much of the data base. 

CHECKR performs consistency checks on the input data. 

This test allows the option of generating a fluid and 
structure portion of the same model at different times. 

Both models may then be accessed from the SAVE file and 
fluid/structure interfaces defined. 

FLSTRl generates the fluid/structure interface data. 

Sequence of subroutines required to generate SHELL models. 
GETLIN and GRDGEN generate grids, QEL0UT elements, and 
STRNGR stringers. 

For the S0LID model, SEL0UT generates elements. 

A general fluid (GFLUID) is defined by GENS0L and the 
associated elements are again found by SEL0UT. 

TFULL (a combined SHELL and S0LID model) is simply the combi- 
nation of the two portions, plus a call to FLSTRl to generate 
fluid/struct are boundaries. 

SHLCAP defines both the grid and element data for a shell 
cap. 

S0LCAP performs grid and element definitions for a solid cap. 

A snell cap containing fluid may be modeled with this routine. 

Two plot routines are available. PL0T generates three- 
dimensional plots, while PL0T2 generates two-dimensional 
pseudo-developed plots for some models. 


Detailed descriptions of these routines and 
their associated utilities are described in 
subsequent sections. 
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2, FILES AND COMMON BLOCKS 


2.1 F0RTRAN FILE DESCRIPTION 

MESHGEN requires the use of six F0RTP^ files in addition to standard 10 units 
5, 6 and 7. The file numbers are found in 

C0MM0N/SYS2/IU1 , IU2 , IU3 , IU4 , IU5 , IU6 

These values are set by Block Data MESHBD and may be changed if desired. 
Operations on each of these files are described below. 

File Usage 

lUl Internal file used to store data passed to plot routines 

IU2 Copy of input stream with MESHLAN data deleted 

IU3 SAVE file for models (usually catalogued by user) 

IU4 Scratch file used in computing fluid/structure interface 

xU5 Available for defining external plot file for SC4020 plotter 

hardware 

IU6 Copy of input stream after conversion to BCD 

The actual formats of the data on each of these files is summarized below. 


2.1.1 File lUl 


Record 


Contents 


1 


2 


3 


ITfPE - tyne of model: 1 = shell, 2 * solid 

NAME - 2 -word BCD M0DEL name 

NR - number of radial stations 

NT - number of circumferential stations 

NZ - number of axial stations 

ICL0SR - flag signifying radial closure 
ICL0ST - flag signifying circumferential closure 

NFLAG 

ANG1,ANG2 - plot view angles 
D0UT 


4 thru n 


GID - point ID 

R,T,Z - grid coordinates (cylindrical) 



2.1.2 Files IU2 -.ad IU6 


The files are simply copies of the input stream, BCD, 80 characters per 
record. 

2.1.3 File IU3 

Record Contents 

1 NAME - 2-word BCD M0DEL name 

2 thru n Bulk Data images - 80 characters per record 

In addition, the boundary element dimensions are saved to find boundaries 
at a later time. The format of such records is: 

Word Contents 

1 BCD String '****' 

2 Element ID 

3,4 ZMIN and 2MAX for element 

5,6 TMIN and TMAX for element 

2.1.4 File IU4 

This file is used internally to find the boundaries when all data is present, 
i.e., SHAPE=TFULL or CAPE. Its format is the same as that described above. 

2.1.5 File IU5 

The format of this file is determined by the SC4020 plot routines. 

2.2 COMMON BLOCKS 

2.2.1 C0MM0N/C0NTRL/M0DNAM (2 ) , SHAPE , B0UND ( 3 ) , NEW 

M0DNAM - Array containing data from M0DEL command 

1,2 M0DEL name, BCD 

3 SAVE flag - 0 ;io SAVE 
1 SAVE 

SHAPE - Type of geometric shape being analyzed from SHAPE command 

= 1 for S.IELL 

- 2 for S0LID 

» 3 for VRUNCATED S0LID 
=« 4 for CONTAINED S0LID 
= 5 for CAPSHELL 

- 6 for CAPS0LID 
= 7 for CAPB0TH 
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, B0UND - Array storing boundary data from B0l)NDAilY command 

i 

• 1 - Type of boundary specification » -1 for INTERP0LATED TABLE 

- 0 for FUNCTI0N 

- 1 for TABLE 

2 •* Boundary ID, integer 

3 - Fluid level for SHAPE-4 or 7 

NEW - Flag to define status of SAVE file. Set to 1 if file is new. 

2.2.2 C0MM0N/ELTAB/ELEM(1O,8) 

ELEM - Array holding element definition for jp to 10 different elements 

per case. Column values are: 
f 

1 - Element type coda 

types 1-8 are quad plates 

types 9-12 are fluids 

types 13-17 are triangular plates 

2 - Initial element ID, integer 
3,4 - Not used 

5 - Element property ID, integer 

6 - Material orientation angle, a, real 

7.8 - Used to denote THICKNESS VARIES command 

2.2.3 C0MM0N/MLINE/NRL,NTL,NZL,R(2OO,9) ,T(200,9) ,Z(200,0) 

This common block is need co store mesh line data generated by expanding 
and merging the zone deflnlti'>n data in /MESH/. 

NRL,NTL,NZL - Number of mesh lines in each coordinate direction, integer 

R,T,Z - Arrays holding expanded mesh line values. Each column of the 
array holds : 

1 - Expanded coordinate values, real 

2 - Expanded coordinate grid ID for shell integer 

3 - Input coordinate system ID for shell in*'eger 

4 - Output coordinate system ID for shell integer 

5 - Not used 

6.7.8 - Same as 2,3,4 for the solid 
9 - Not used 

2.2.4 C0MM0N/MESH/G1,C2,0PR0P,X1(5,2O) ,X2(5,20) ,X3(5,20) ,NPT(3) ,STRING(6) , 

!. SL0C(2,25) 

G1 - Initial grid ID for shell, lnteg.^r 

G2 - Initial grid ID for solid, integer 
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- Overall property ID, integer. 


X1,X2,X3 - Arrays storing zone definition values for each coordinate 

direction. The first subscript refers to the number of the 
STEP or DIVIDE command in the ibh direction, l.e., a maximum 
of 5 STEP or DIVIDE 's may appear in each direction. The 
words within each entry are: 

1. First value of coordinate on STEP command, real. 

2. Second value of coordinate on STEP command, real. 

3. Number of elements in the zone from STEP or DIVIDE 
command, integer. 

4. Shell flag - set to 1 if plate elements are defined for 
the zone, integer. 

5. Input coordinate system ID for the shell zone, INSYS, 
integer. 

6. Output coordinate system ID for the shell zone, 0UTSYS, 
integer. 

7. Grid point ID Increment for shell zone from NUMBER command, 
integer. 

8. Element ID increment for shell zone from NUMBER command, 
integer. 

9. Pointer into ELEM array of /ELTAB/ defining the type of 
elements in the shell zone, integer. 

10. Number of element types defined in the shell zone, 
integer. 

11. Shell zone property ID from ZPR0P command, integer. 

12. Mrterial orientation angle, a, for shell zone from ZPR0P 
command , real . 

13-19. Same as 4-10 for the solid zone. 

20. Solid zone property ID from the ZPR0P command, integer. 
NPT(I) - Nimber of zones defined in each coordinate direction. 

STRING - 1 - Type of stringer, 0 for BAR, 1 for R0D. 

2 - Stringer property ID, integer. 

3 - Initial stringer ID number if 0 direction, integer. 

4 - Niunber of stringer stations in Z direction from AL0NG 

command, integer. 

5 - Number of stringer stations in 0 direction from AL0NG 

command, integer. 

6 - Initial stringer ID number if Z direction, integer. 

SL0C - Array containing stringer stations from AL0^’G command <n the Z 
and 0 directions. 



3. UTILITY SUBROUTINES 


3.1 MAPFNS (MACHINE WORD FUNCTIONS) 

3.1.1 Entry Points ; LSHIFT, RSHIFT, ANDF, 0RF, X0RF, C0MPLF 

3.1.2 Purpose 

To perform basic computer word manipulations by standard binary digit (bit) 
operations. The manipulations ata performed over the complete memory word 
length for the particular hardware. 

3.1.3 Calling Sequence 

All machine word functions are executed as F0RTRAN integer function sub- 
routines with integer arguments. 

3.1.4 Method 


The method employed within each function will be described following the 
separate function examples. 

3.1.5 Entries 


K = LSHIFT (I,N) 

The entire bit structure of word I is shifted left N places and the resulting 
word replaces word K. Word I is unchanged. High-order bits shifted out are 
lost. Zeros are supplied to vacated low-order positions. The shift is logical; 
no special provision is made for the sign position. 

K = RSHIFT (I,N) 

The entire bit structure of word I is shifted right N places and the resulting 
word replaces word K. Word I is unchanged. Low-order bits shifted out are 
lost. Zeros are supplied to vacated high-order positions. The shift is logical 
no special provision is made for the sign position. 

K = ANDF (I,J) 

A logical product of the bits within word I and word J is formed and stored 
into word K. Words I and J are unchanged. 

K “ 0RF (I,J) 

A logical sum of the bits within word I and word J is formed and stored into 
word K. Wordj I and J are unchanged. 
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K = X0RF (I»J) 

The modulo-two sum (exclusive or) of the bits within word I and word J is 
formed and stored into word K. Words I and J are unchanged. 

K = C0MPLF (I) 

The ones complement of the bits within word I is formed and stored into 
word K. Word I is unchanged. 

3.1.6 Design Requirements 

NAPFNS is written in assembly language. 


\ 
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3.2 XRCARD (FREE-FIELD CARD DATA CONVERSION ROUTINE) 

3.2.1 Encvy Point ; XRCARD 

3.2.2 Purpose 

To Interpret free-field card input data as follows: 

1. Identify BCD alpha and numeric data fields as they are converted 
and placed in the user's buffer. 

2. Flag and output special data field delimiters. 

3. Convert BCD numeric fields to binary Integer or binary floating 
point. 

4. Indicate when the data extends beyond one 72 column card. 

3.2.3 Calling Sequence 

CALL XRCARD (0UTBUF,L,INBUF) 

Where: 

0UTBUF = The buffer which is to contain the converted card image. 

L = The length of 0UTBUF available to XRCARD. 

INBUF = The buffer containing the card image to be converted. 

3.2.4 Method 

XRCARD' s design is based on the necessity of having to function on a variety 
of computing machines having a variety of computer word structures, and a 
variety of differences in hollerlth handling imposed by differing F0RTRAN 
compilers . 

XRCARD analyzes the twenty hollerith words input thro< INBUF as follows: 
Data Field Delimiters 


t* 

•* tr 


Type A: The following symbols signify the end of an alpha field or numeric 

field on the card. As these symbols are encountered, they will be 
flagged and placed in the output buffer to aid the user in Identify- 
inr the data. 

( LEFT PAREN 
/ SLASH 
= EQUAL 
* ASTERISK 


■« 

* 1 


Type B: The following symbols are identical to those listed above except 

that the symbol not flagged or placed in the output buffer: 
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, COMMA 
) RIGHT PAREN 

When successive type A or type B delimiters are encountered, a null 
field indication (two BCD blank words) is output. A null field is 
generated for each successive delimiter. A null field is also gener- 
ated when a type A or type B delimiter is followed by a $ indicating 
the end of data condition. 

Type C: The following symbol is identical to the COMMA except that no null 

field indication is output when they are encountered in succession. 

BLANK 


End of Data Indication 

There are three means by which end-of-data may be specified on the card: 

• The last data field ends in column 72, or is followed by blanks out 
through column 72; 

• $ is encountered, after which comments may be included out to column 
80; or 

• Continuation cards ending in (, /, = or , will result in a continuation 
flag (0 mode word)). 


Format of Output Data 

A mode word, N, is placed in the output buffer to distinguish between BCD 
data and numeric data. 


Numeric Mode Word : A new mode word is output each time a numeric field 

is converted and output. (All numeric mode words are negative.) 

N « -1 Integer data (1 data word) 

= -2 floating point single precision (1 data word) 

= -4 floating point double precision (2 data words) 

N indicates the type of numeric data and where to look for the next mode 
word. 


Alpha Mode Word : When processing alpha data, only one mode word is output 

for successive alpha fields, i.e. , an alpha mode word will never follow 
another alpha mode word. 

N = The number of successive alpha fields encountered on the card. 
Each alpha field consists of two 4-character computer words 
(left adjusted). Thus, N can be used to compute the location 
of the next mode word. 

The type A delimiters are output as alpha data and are 'covered' by the 
alpha mode word. Since data output in the alpha mode must consist of two 
words, a type A delimiter will appear as: 

Word 1 = Delimiter flag, all bits of the word are on. 

Word 2 = BCD delimiter, left adjusted, followed by BCD blanks. 
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End-of-Data ; The end-of-data flag is placed last in the output buffer 
and appears in place of an expected mode word. There are two end-of-data 
flags: 

• A word with all bits off, indicating that more data is to follow 
on a continuation card. 

• A word will all bits on except for the sign, indicating that no 
more data is to follow for this card type. 

3.2.5 Design Requirements 

An alpha field must be eight characters or less. Long alpha fields will be 
truncated to eight characters. 

All data must be placed in card columes 1-72 . 

A data field may not be split between two cards. 

The specification of all numeric data fields must conform to F0RTRAN IV 
standards. 


3.3 RE2AL (REAL NUMBER TO ALPHANUMERIC) 


3.3.1 Entry Point ; RE2AL 

3.3.2 Purpose 

To convert a single precision number to its BCD string representation. 

3.3.3 Calling Sequence 
CALL RE2AL(RE,ALPH) 

where : 

RE - single precision real number - real - input. 

ALPH - BCD output (2 words) 

3.3.4 Method 

The output value is the floating or exponential form depending on which gives 
greater significant digits. Round off is based on the number of significant 
digits being output. 



3- A INT2AL 


3.4.1 Entry Point : INT2AL 

3.4.2 Purpose 

To convert an integer value to its BCD string representation. 

3.4.3 Calling Sequence 
CALL INT2AL(NUMB,ALPH) 

where ; 

NUMB - number to be converted to BCD string - integer - input 
ALPH - two-word arriiy containing BCD string - BCD - output 

3.4.4 fiethod 

Each digit of NUMB is isolated and the BCD string is built by appropriate 
shift and logical "OR" operations. If the integer value is greater than 8 
digits, the string returned in ALPH is all "question marks" (?). 


5-14 


#■ 



3.5 PAGE (PAGE HEADING) 


3.5.1 Entry Point ; PAGE 

3.5.2 Purpose 

To provide a standard page head'*ng foi MESHGEN output. 

3.5.3 Calling Sequence 
CALL PAGE 

C0MM0N/SYSTEM/NLPP ,NLINE , IPAGE 

NLPP - Maximum number of lines per page - integer. 

NLINE - Number of data lines on previous page - LINE is set to 
zero by PAGE. 

IPAGE - Current page number - Increased by 1 on each call to PAGE. 
C0MM0N/0UTPUT/NTITL(32) ,ITITL1(32) ,ITITL2(32) ,ITITL3(32) 

3.5.4 Method 

PAGE writes a standard 4 line heading from NTITL, ITITLl, IT1TL2, ITITL3. 



3.6 MESAGE (MESSAGE WRITER) 


3.6.1 Entry Point ; MESAGE 

3.6.2 Purpose 

To queue nonfatal messages during the execution of a module; and for fatal 
messages give a core dump (CALL PDUMF) , print the message queue (CALL MSGWRT) , 
and call PEXIT. 

3.6.3 Calling Sequence 
CALL MESAGE (IC0DE,L,M,N) 

where : 

IC0DE - Internal message number - integer - input. 

L,M,N - Values to be Included in the error message - mixed - input. 

3.6.4 Method 


Both fatal and nonfatal error messages are printed as they occur. If a 
fatal error occurs, the lERR flag is set and return is made to the calling 
routine where appropriate action must be taken. 



3.7 BUG 


3.7.1 Entry Point : BUG 

3.7.2 Purpose 

To allow for diagnostic output during program development. 


3.7.3 Calling Sequence ; 

CALL BUG(IDENT,IDNUM, ARRAY, NWDS) 
where : 

IDENT - Four-character Identifier to label output - BCD - input. 

IDNUM - Number to identify output - integer - input. 

ARRAY - Starting address of values to be output - input. 

NWDS - Number of words to be output - Integer - input. 

C0MM0N/DIAG/IDIAG 

3.7.4 Method 

If the diagnostic flag IDIAG Is on, the routine will print the specified 
data. The type of data and output formats are determined by the subroutine. 
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3.8 APRXEQ (APPROXIMATELY EQUAL FUNCTION) 


3.8.1 Entry Point ; APRXEQ 

3.8.2 Purpose 

To determine whether two single-precision real values are approximately 
equal within a given tolerance. 

3.8.3 Calling Sequence 

K = APRXEQ(A,B) [LOGICAL FUNCTION] 
where : 

AyB - Numbers to be compared for approximate equality - real - Input. 
C0MM0N/SYS1/T0LER 
where : 

T0LER - Approximate equality tolerance - real. 

3.8.4 Method 

The following algorithm Is used to determine approximate equalxty: 

; f A=B; then APRXEQ* . TRUE . 
else 

U (Aj«OA|~l<T0LER)V 

(B}«0A|^| <T0LER) ; then APRXEQ-. TRUE, 
else 


APRXEQ-. TRL»E. 



3.9 DECIDE (D0F DECODER) 


3.9.1 Entry Point ; DEC0DE 

3.9.2 Purpose 

To translate the normal NASTRAN degree-of~freedom code (string of integers l-’6 
with no imbedded blanks) to a bit representation. 

3*9*3 Calling Sequence 

K - DEC0DE(D0F) (INTEGER FUNCTION) 
where: 

D0F - NASTRAN degree-of-freedom code - Integer - input. 

3. 9. A Method 

The digits of D0F are Isolated and then a word Is constructed where the 1 
bit is on if the digit 1 appears in D0F. For example: 

D0F « 126; DEC0DE « 100011 
D0F » 123A56; DEC0DE - 111111 
D0F - 256; DEC0DE = 110010 
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3.10 ENC0DE (D0F ENCODES) 


3. -10.1 Entry Point ; lOTEGER FUNCTI0N ENC0DE 

3.10.2 Purpose 

To translate a bit pattern to the NASTRAN degree-of •'freedom code. 

3.10.3 Calling Sequence 

K - ENC0DE(D0F) (INTEGEP FACTION) 

where : 

D0F - Bit string to be coded - integer - Input. 

3.10.4 Metho d; 

Bits are extracted one at a time and the D0F string Is generated, 
of this process are: 

011011 -> 1245 
111111 “> 123456 
000011 12 


Examples 


This routine Is the inverse of DEC0DE. 



3.11 INTERP 


3.11.1 Entry Point ; INTERP 

3.11.2 Purpose 

To perform a linear Interpolation on a tabular function r = f(z). 

3.11.3 Calling Sequence 

R = INTERP (Z) (REAL FUNCTION) 
where 

Z - value of axial coordinate at which radial value is desired - real - 
inpit. 

C0MM0N/BXAB/ 


3.11.4 Method 


A linear interpolation is performed to find the radial value. An error 

condition exists if Z < Z . or Z > Z 

min max 
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3.12 GETVAL 


3.12.1 Entry Point ; GETVAL 

3.12.2 Purpose 

To return a value of radial coordinate for a given axial value. 

3.12.3 Calling Sequence 

R = GETVAL(Z) [Real Function] 
where : 

2 - axial value ; inpi t , real 
C0MM0N/C0NTRL/ 

3.12.4 Method 

A test is made to determine whether the boundary has been specified by a 
FUNCTI0N or a TABLE. The appropriate call is then made to either FUNC or 
INTER? . 
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3.13 FUNC 


3.13.1 Entry Point ; FUNC 

3.12.2 Purpose 

Evaluates the incomplete quadratic function 

a^z^ + 352 + agr^ P a^r - + a^ 

for r given a particular value of z. 

3.12.3 Calling Sequence 

R = FUNC(Z) (Real Function) 
where : 

Z - value of axial coordinate at whicl radial coordinate is desired 
input , real . 

C0MM0N/SHP/A(9) 

A(I) - array containing quadratic coefficients; real. 

3.13.4 Method 

The following algorithm is used to find the value of r: 

If a, - a_ = 0 then ERR0R 

o / 2 2 

If a^ = 0 then FUNC = (ag + a^ - a^z - a^z)!a^ 

2 2 2 
else DISC = - 4ag(a^z + a^z ~ 

if DISC 0.0 then ERR0R 

R1 = (-a^ + v^ISC)/2ag 
R2 = (-a^ - t'^I^)/2ag 
FUNC = max(Rl,R2) 
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3.14 DIVLIN 


3.14.1 Entry Point ; DIVLIN 

3.14.2 Purpose 

To divide a boundary curve into n approximately equal segments. 

3.14.3 Calling Sequence 
CALL DIVLIN (NPT.ZP) 

where : 

NPT - the number of equally-spaced points - input - integer. 

ZP - array of values at equally-spaced points (axial locations) 
output - real. 

C0MM0N/BTAB/ C0MM0N/SHP 
C0MM0N/C0NTRL/ 


3.14.4 Method 

The SHAPE flag in the /C0NTRL/ common blocK is checked to determine whether 
the boundary is defined by a function or table. In either event the length 
ot the curve is approximated by the summation of line segments defining the 
curve. This length is then divided into NPT equal segments and interpolation 
is performed to give the axial values of these points. 


5- 24 


3.15 PIFl 


3.15.1 Entry Point ; PIFl 

3.15.2 Purpose 

To perform a linear interpolation of the curve length function. 

3.15.3 Calling Sequence 
CALL PIF1(Z,L,N.LP,2P) 

where : 

Z - table of axial values - real - input. 

L - table of arc length at each Z - real - input. 

N - number of points in tables - integer - input. 

LP - particular value of length for which the axial value is 

desired - real - input. 

ZP - Interpolated value of Z at LP - real - output 
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4, DATA GENERATION AND 10 SUBROUTINES 


4.1 FRE0UT 


^•1*1 Entry Point ; FRE0UT 

4.1.2 Purpose 


To generate fluid-free surface CFFREE Bulk Data cards. 

4.1.3 Calling Sequence 
CALL FRE0UT (ID.IFACE) 

where: 

ID - Fluid element ID number - integer - input. 

IFACE - Face identification for free surface - integer - input. 

C0MM0N/L0ADS/ 

4.1.4 Method 


Subroutine INT2AL is called for each of the input values and the Bulk Data 
card image is generated. The image is then written to the output file. In 
addition, if the SAVE or PUNCH flags are on, the a.mage will be written to 
the SAVE file and/or punched. 
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4.2 GRD0UT 

4.2.1 Entry Point ; GRD0UT 

4.2.2 Purpose 

To generate GRID Bulk Data cards. 

4.2.3 Calling Sequence 

CALL GRD0UT(ID,IN,R,T.Z,0UT,SPC) 
where : 

ID - grid point ID number - input - Integer. 

IN - input coordinate system ID - input - integer. 

R,T,Z - grid point coordinates - input - real. 

0UT - output coordinate system ID - input - integer. 

SPC - permanent single point constraint code - input - integer. 


4.2.4 Method 


Subroutine INT2AL is called for each input value and the Bulk Data card image 
is generated. The value SPC must be encoded back to NASTRAN form. In addition, 
if the PL0T flag is set, the grid ID and coordinates are written to the plot 
file. If the SAVE or PUNCH flags are on, appropriate action is taken. 


< 


I 
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4.3 HEX0UT 


4.3.1 Entry Point ; HEX0UT 

4.3.2 Purpose 

To generate CFHEXl and CFHEX2 Bulk Data cards. 

4.3.3 Calling Seque nce 

CALL HEX0UT(ITYPE,E,P ,G1 ,G2 ,G3 ,G4 ,G5 ,G6 ,G7 ,G8) 
where : 

ITYPE - flag to determine type of card - input - integer. 
= 9 CFHEXl 
= 10 CFHEX2 

E - element ID number - input - integer. 

P - property ID number - input - integer. 

G1-G8 - grid point ID number - input - integer. 


4.3.4 Method 

Subroutine INT2AL is called for each of thu input values and the Bulk Data 
card image is generated. The image is then written to the output file. In 
addition, if the SAVE or PUNCH flags are on, the image will be written to the 
SAVE file and/or punched. 
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4.4 PL0TEL 


4.4.1 Entry Point ; PL0TEL 

4.4.2 Purpose 

Generates PL0TEL Bulk Data cards for the free surface of a fluid model. 


4.4.3 Calling Sequence 

CALL PL0TEL(ID1,G1,G2,ID2,G3,G4) 
where ; 

ID1,ID2 - PL0TEL element ID numbers - input - integer. 
G1,G2,G3,G4 - grid point ID's defining the PL0TEL - input - integer. 


4.4.4 Method 

Subroutine INT2AL is called for each of the input values and the Bulk Data 
card image is generated. The image is then written to the output file. In 
addition, if the SAVE or PUNCH flags are on, the image will be written to 
the SAVE file and/or punched. 


4.5 PLTHED 


4.5.1 Entry Point ; PLTHED 

4.5.2 Purpose 

To write three header records on the scratch file used by the plotting 
routines, 

4.5.3 Calling Sequence 
CALL PLTHED(ITYPE,IP0S) 

where : 

ITYPE - type of plot that will be generated 
= 1 for shell 
= 2 for solid 

IP0S - code for positioning of file 
= 1 rewind 
=» 2 no rewind 

4.5.4 Method 


Twelve words are written as three four-word header records to the plot file. 
(See Section 2.1 for the format of this file.) 



4.6 PNCHER 


4.6.1 Entry Point ; PNCHER 

4.6.2 Purpose 

To punch the Bulk Data for a model that resides on the SAVE file. 

4.6.3 Calling Sequence 
CALL PNCHER 

C0MM0N/I0DATA/JUNK(5) ,NAME(2) 
where : 

NAME - 2-word array containing the model name on the SAVE file. 

4.6.4 Method 

The SAVE file is accessed and a search Is made for model NAME. If NAME is 
found, all Bulk Data Images are punched. If NAME does not exist, the PUNCH 
command is ignored. 


4.7 STR0UT 


4.7.1 Entry Point ; STR0UT 

4.7.2 Purpose 

To generate CH0D and CBAR Bulk Data cards. 


4.7.3 Calling Sequence 

CALL STR0UT(IC,EID,PID,G1,G2) 
where : 

IC - code for element type - input - Integer 

- 1 for CBAR 

- 2 for CR0D 

EID - element ID number - input - integer 
PID - property ID number - input - Integer 
G1,G2 - grid point ID numbers - input - integer 


4.7.4 Method 


Subroutine INT2AL is called for each of the input values and the Bulk Data 
card image is generated. The image is then written to the output file. In 
addition, if the SAVE or PUNCH flags are on, the image will be written to 
the SAVE file and/or punched. 



A. 8 TET0UT 


4.8.1 Entry Point ; TET0UT 

4.8.2 Purpose 

To generate. CFTETRA Bulk Data cards. 

4.8.3 Calling Sequence 

CALL TET0UT(E,P,G1,G2,C3,G4) 
where : 

E - element ID nom>'°'' - input - Integer. 

P - property IP number - ^ out - Integer. 

G1-G4 - grid point ID numbers - input - integer. 

4.8.4 Method 

Subroutine INT2AL is called for e?v i of the input values and the Bulk Data 
card image is generate.'. The image is then writte*' to the output file. In 
addition, if the SAVE or PITNCH flags are on, the image will be written to 
the SAVE file and/or punched. 
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4.9 WED0UT 


4.9.1 En try Point ; WED0UT 

4.9.2 Purpose 

To generate CFWEDGE Bulk Data cards. 

4.9.3 Calling Sequence 

CALL WED0UT(E,P,G1,G2,G3,G4,<'5,G6) 
where ; 

E - element ID number - input - integer. 

P - propel tj- ■ 3 number - input - integer. 
G1-G6 - grid point ID numbers - input - integer. 


4.9.4 Method 

Subroutine INT2AL is called for each of the input values and the Bulk Data 
card image is generated. The image is then written to the ou'put file. In 
addition, if the SAVE or PUNCH flcgs .•’re on, the image will be written to 
the SAVE file and/or punched. 
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4.10 WRTHED 

4.10.1 Entry Point ; WRTHED 

4.10.2 Purpose 

To write a header record on the SAVE file if the SAVE option is elected for 
a given model. 

4.10.3 Calling Sequence 
CALL WRTHED 
CC;.L.0N/C0NTRL/ 


* 
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4.11 CCARD 


4.11.1 Entry Point ; CCARD 

4.11.2 Purpose 

To generate the Bulk Data card continuation field. 

4.11.3 Calling Sequence 

S 

CALL CCARD (NUM.0UT, PREFIX) 
where : 

NUM - numeric portion of continuation field - input - Integer. 

0UT - 2-word array containing continuation field - output - BCD. 
PREFIX - 2-character prefix for continuation field - input - BCD. 

4.11.4 Method 

Logical shifting and bitwise logical "OR" functions are performed to pack 
the prefix onto the BCD representation of NUM. 
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4.12 PREPAS 


4.12.1 Entry Point: PREPAS 

I 

4.12.2 Purpose 

To convert the input data stream from EBCDIC to BCD. 

4.12.3 Calling Sequence 
CALL PREPAS 

, 4.12.4 Method 


PREPAS performs a character-by-character conversion of the input stream and 
writes the input to unit 1U6. 

4.12.5 Design Requirements 

PREPAS is an IBM-dependent routine and must be altered to simply copy the 
input data on other systems. 
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4.13 DATAEX 


4.13.1 E ntry Point ; DATAEX 
4 •13. 2 Purpose 

To extract the MESHGEN functional data from the input stream in a pre-pass 
operation before mesh generation. 

4.13.3 Calling Sequence 
CALL DATAEX 
C0MM0N/B0UND/ 

4.13.4 Method 

The input file is read until the $DATA card is encountered. Each FUNCTI0N 
or TABLE definition is extracted and loaded into /B0UND/. When a hardware 
end-of-data is sensed on the input unit , a REWIND is performed and the 
MESHLAN sequence is copied to a scratch unit and all cards after and Including 
$DATA are eliminated. 
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4.14 INIT 


4.14.1 Entry Point ; INIT 

4.14.2 Purpose 

To extract the boundary function or table data for the current case from the 
tables of all data input for the current job. 


4.14.3 Ca lling Sequence 
CALL INIT 

C0MM0N/B0UND/FT(1O,1O) ,TT(500) 

C0MM0N/SHP/A(9) 

C0. M0N/BTAB/NPT,TABL(5OO) 

/B0UND/ contains tables of all input functions and boundary tables 
/SHP/ current function data 
/BTAB/ current table data 

4.14.4 Method 

This routine extracts the current boundary data from the tables FT and TT. 
It then vrites a summary of the data and copies it into A or TABL. 
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4. IS ZER^UT 


4.15.1 Entry Point ; ZER0UT 


4.15.2 Purpose 

To initialize the basic common blocks before executing each case in a 
MESHGEN run. 


4.15.3 Calling Sequence 


CALL ZER0UT 

C0MM0N/MESH/ 

C0MM0N/MLINE/ 

C0MM0N/C0NTRL/ 

C0MM0N/ELTAB/ 


C0MM0N/ SYSTEM/ 
C0MM0N/I0DATA/ 
C0MM0N/PSPC/ 
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4.16 BLOCK DATA MESHBD 


4.16.1 Entry Point ! MESHBD 

4.16.2 Purpose 

To initialize machine-dependent parameters held in common blocks /SYSl/ 
and /SYS2/. 
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5. LEXICAL ANALYSIS SUBROUTINES 


5.1 LEXER 

5.1.1 Entry Point ; LEXER 

5.1.2 Purpose 

To perform the lexical analysis of MESHLAN statements* and to Initialize 
control arrays for each M0DEL. 

5.1.3 Calling Sequence 
CALL LEXER 


5.1.4 Method 


The logical flow of LEXER Is shown In Figure 2. The routine determines the 
level of the MESH6EN command and calls the processing routine. 
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Call ZER0UT 



5 

No 

f 

Get Keyword Token 




ENDJ0B- .TRUE. 

> 



C^TURN^ 



Error: Invalid Keyword 


FIGURE 2. LOGICAL FLOW OF LEXER 
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5.2 LEVELl 


5.2.1 Entry Point ! LEVELl 

5.2.2 Purpose 

To translate and process the primary MESHLAN commands: 

M0DEL, GE0METRY, MESH, GRAV, PL0T, PLT2, PUNCH, ENDM, DTAG, 

FIND, MACH, and PLTH*. 

5.2.3 Calling Sequence 
CALL LEVELl (P) 

where : 

P - Primary command ID - Integer - Input. 

C0MM0N/MESH/ C0MM0N/MACH0N/ 

C0MM0N/L0ADS/ C0MM0N/I0DATA/ 

C0MM0N/BFLSTR/ C0MM0N/C0NTRL/ 

5.2.4 Method 

The primary commands are translated and values Inserted Into the appropriate 
common blocks. 


*PLTHEAD is a facility-dependent command for MSFC. 
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5.3 LEVEL2 


5.3.1 Entry Point : LEVEL2 

5.3.2 Purpose 

To translate and process the secondary MESHLAN commands: 

SHAPE, B0UNDARY, STEP, DIVIDE, and 0PR0P. 

5.3.3 C alling Sequence 
CALL LEVEL2(P,S) 

where : 

P,S - Primary and secondary command ID's - Input- Integer. 

C0MM0N/C0NTRL/ 

C0MM0N/MESH/ 

5.3.4 Method 


Tests are performed to check the context of the secondary commands. The 
commands are then translated and' the values Inserted Into the common blocks. 



5. A LEVEL3 


5.A.1 Entry Point ; LEVEL3 

5.4.2 Purpose 

To translate and process tertiary MESHLAN commands: 
SHELL, S0LID, ZPR0P, and STRINGER. 


5.4.3 Calling Sequence 
CALL LEVEL3(P,S,T) 
where : 

P,S,T - Primary, secondary and tertiary command ID’s - Integer - input. 
C0.'^LM0N/MESH/ 


5.4.4 Method 

Tests are performed to check the context of the tertiary commands. After 
translation, the appropriate values are inserted into /MESH/. 


5“ 46 



5.5 LEVELA 


5.5.1 Entry Point ; LEVELA 

5.5.2 Purpose 

To translate and process the quaternary MESKLAN coinmands : 

NUMBER, INSYS, 0UTSYS, FIX, ELEMENTS, and AL0NG. 

5.5.3 Calling Sequence 
CALL LEVELA (P,S,T,Q) 

where : 

P,S,T,Q ~ primary, secondary, etc. command ID's - input - integer. 
C0MM0N/MESH/ C0MM0N/ELTAB/ C0MM0N/PSPC/ 

5. 5. A Method 

Tests are performed to check the syntax of the LEVELA commands. The values 
extracted are then Inserted into the appropriate common block. 
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5.6 LEVELS 


5.6.1 Entry Point : LEVELS 

5.6.2 Purpose 

To process the quintary MESHLAN commands: PR0PERTY and THICKNESS. 

5.6.3 Calling Sequence 
CALL LEVEL5(P,S,T,Q1,Q2) 

where : 

P,S,T,Q1,Q2 - command ID for primary, secondary, etc. command sequence. 
C0MM0N/ ELTAB / ELEM ( 10 , 8 ) 

5.6.4 Method 

Tests are performed to check the syntax of the LEVELS commands. Extractei 
values are then placed in the /ELTAB/ common block. 



5.7 CHECKR 


5.7.1 Entry Point : CHECKR 

5.7.2 Purpose 

To test for the allowability and consistency of the processed MESHLAN 
sequence . 

5.7.3 Calling Sequence 
CALL CHECKR 

C0MM0N/MESH/ C0MM0N/ELTAB/ 

C0MM0N/L0ADS/ C0MM0N/C0NTRL/ 

C0MM0N/BFLSTR/ 
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5.8 GETT0K 


5.8.1 Entry Point ; GETT0K 

5.8.2 Purpose 

To return the next token of a MESHLAN command and its type. 

5.8.3 Calling Sequence 
CALL GETT0K 

COMMON/LEXrC/TC0DE ,T0K(2) 
where : 

TC0DE - Is the token type (see Section 3.2) 

T0K - is the actual value of the token 


5-50 



5.9 GETCRD 


5.9.1 Entry Point ; GETCRD 

5.9.2 Purpose 

To process the MESHLAN commands for lexical analysis. 

5.9.3 Calling Sequence 
CALL GETCRD 

C0MM0N/CARD/T0KEN(3,1OO) ,NT0K 
where : 

T0K - Array containing the types and values of each token extracted 
from the Input stream. The contents are: 

Row Contents 

1 Token type code 

1 - BCD 

2 - Integer 

3 - Real 

4 - Delimiter 

5 - End-of-card 

2 Value of token 

5.9.4 Method 

When GETCRD Is called. It reads the next card image In the input stream. 
XRCARD is then called to extract the tokens from the card. This routine 
then processes the XRCARD output to a simpler form to be used by the lexical 
routines. Although this function Is purely overhead. It Is of great use In 
simplifying coding and understanding of the lexical routines LEVELl - LEVEL5 . 
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6 . 


PLOTTING subroutines 


6.1 PL0T2 

6.1.1 Entry Points ; PL0T, PL0T2 

6.1.2 Purpose 

To drive the various plotting subroutines. 

6.1.3 Calling Sequence 

CALL PL0T2 (FRAME) - 2-D plots 

CALL PL0T (FRAME) - 3-D plots 

where : 

FRAME - last frame number, zero initially 


6.1.4 Method 

The two entry points set a flag, IPLT23, that determines if 2-D or 3-D plots 
are to be made. The title is printed, then PLTRD is called to initialize the 
on-core arrays. If the plot can be made, various subroutines are called to 
do each aspect. 
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6.2 PLTNFR 


6.2.1 Entry Point ; PLTNFR 

6.2.2 Purpose 

For SC-4020 plots, skip to a new frame, write the plot titles, and increment 
the frame number. 

6.2.3 Calling Sequence 
CALL PLTNFR (FRAME) 

where : 

FRAME - input/output frame number 

6.2.4 Method 

If the frame number is £ 0, the SC-4020 is initialized. The plot titles 
are generated from data in /TITLE/ and /PLTDAT/ . 

6.2.5 Additional Subroutines 

SC-4020 routines LABLV, CAMAY, FRAMEV, and PRINTV are called. 
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6.3 PLTRD 


6.3.1 Entry Point ; PLTRD 

6.3.2 Purpose 

For SC-4020 plots, create the in-core arrays of grid data. The data is 
rotated as per user request and converted to rasters. Each type of element 
input causes a resorting and/or addition of points on the input file. 


6.3.3 Calling Sequence 

CALL PLTRD(IFIL,IGRL,IGR,RGR) 


where : 


IFIL - input file number, binary records of 4 words. 


IGRL - input, number of words of open core. 


1GR\ 

RGRj 


output array of grid data, 3 words per point 

word 1 = grid ID 
word 2 = X raster 
word 3 = Y raster 


6.3.4 Method 


Each of the element types, defined in the header record of the input file, 
causes a special sort of the grid data output. For closure in Z or theta, 
extra points are introduced. 



6.4 TPLAB 


6.4.1 Entry Point ; TPLAB 

6.4.2 Purpose 

For SC-4020 plots > places the grid identification number to the right of the 
grid point. 

6.4.3 Calling Sequence 
CALL TPLAB (IGR) 

where : 

IGR - 3-word array of grid data 

Word Description 

1 grid ID number 

2 X raster 

3 Y raster 

6.4.4 Method 

The grid set is converted to alphanumeric and printed 32 rasters to the 
right of the point. 

6.4.5 Additional Subroutines 


SC-4020 routine PRINTV is required. 
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6.5 


TPLIN 


6.5.1 Entry Point ; TPUN 

6.5.2 Purpose 

For SC-4020 plots, to call one of the four line-drawing routines. 

6.5.3 Calling Sequence 
CALL TPLIN (IGR) 

where : 

IGR - input 3-word grid ID and raster array 
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6.6 TYPLl 


6.6.1 Entry Point ; TYPLl 

6.6.2 Purpose 

For SC-4020 plocs, to draw the elements tor a shell of revolution model. 

6.6.3 Calling Sequence 
CALL TYPLl (IGR) 

where : 

IGR - 3-word grid ID and raster array 

6.6.4 Method 

The pre-deflned array of IGR set up by PLTRD is scanned to extract raster 
locations. 

6.6.5 Additional Subroutines 
SC-4020 routine LINEV is required. 
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6.7 TYPL2 


6.7.1 Entry Point ! TYPL2 

6.7.2 Purpose 

For SC-4020 plots, to draw the elements for a solid of revolution. 

6.7.3 Calling Sequence 
CALL TYPL2(IGR) 

where: 

IGR - 3-word grid ID and raster array 

6.7.4 Method 

The pre-deflned array, IGR, set up by PLTRD is scanned to extract raster 
locations. 

6.7.5 Additional Subroutines 
SC-4020 routine LINEV is required. 
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6.8 TYPL3 


6.8.1 Entry Point ; TYPL3 

6.8.2 Purpose 

For SC-4020 plots, to draw the elements for a spherical shell cap. 

6.8.3 Calling Sequence 
CALL TYPL3(IGR) 

where : 

IGR - 3-word grid ID and raster array 

6.8.4 Method 

The grid-defined array, IGR, set up by PLTRD is scanned to extract raster 
locations. 

6.8.5 Additional Subroutines 
SC-4020 routine LTNEV is required. 
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6.9 TYPL4 


6.9.1 Entry Point ; TYPL4 
6.9.J? Purpose 

For SC-4020 plots, to draw the elements for a spherical solid cap. 

6.9.3 Calling Sequence 
CALL TYPL4(IGR) 

where ? 

IGR - 3-word grid ID and raster array 

6.9.4 Method 

The pre-deflned array, IGR, set up by PLTRD is called to extract raster 
locations. 

6.9.5 Additional Subroutines 
SC-4020 routine LINEV is required. 
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7. 


COMPUTATIONAL SUa*iOUTINLS 


7 . 1 GETLIN 


7.1.1 Entry Point : GETLIN 

7.1.2 Purpose 

To generate the mesh lines in each coordinate direction and apply permanent 
single-point constraints. 

7.1.3 Calling Sequence 

CALL GETLIN 

C0MM01,, C0NTRL/ 

C0MM0N/MLINE/ 

C0MM0N/MESH/ 

7.1.4 Method 


C0MM0N/PSPC/ 

C0MM0N/SHP/ 

C0MM0N/BTAB/ 


The mesh definition data is accessed from /MESH/. Tuis data is then expanded 
to define all the coordinate values for the model. These are then stored in 
/MLINE/ (see Section 2.2.3). In addition to this expansion, SPCs are also 
applied to each grid line. 


SM 



7.2 QEL0UT 


7.2.1 Entry Point : QEL0UT 

7.2.2 Purpose 

To generate element connection for shell models. 

7.2.3 Calling Sequence 
CALL QEL0UT 

C0MM0N/C0NTRL/ C0MM0N/ML1NE/ 

C0MM0N/MESH/ C0MM0N/ELTAB/ 
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7.3 SEL0UT 


7.3.1 Entry Point ; SEL0UT 

7.3.2 Purpose 

To generate element data for solids of revolution (S0LID, TFULL, GFLUID) 

7.3.3 Calling Sequence 
CALL SEL0UT 

C0MM0N/MESK/ C0MM0N/C0NTRL/ 

C0MM0N/MLINE/ C0MM0N/ELTAB/ 

7.3.4 Method 

The grid point ID numbers and increments are accessed and the appropriate 
solid element connection Bulk Data is generated. 1 1 addition, PL0TEL data 
is generated for free surfaces and element ranges computed for CFLSTR Bulk 
Data. 
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7.4 GENS0L 


7.4.1 Entry Point : GENS0L 

7.4.2 Purpose 

To generate grid points for the generalized solid model (GFLUID) . 

7.4.3 C alling Sequence 
CALL GENS0L 

C0MM0N/C0NTRL/ C0MM0N/MLINE/ 

C0MM0N/BTAB/ C0MM0N/SHP/ 

C0MM0N/MESH/ 

7.4.4 Method 

GENS0L must compute all coordinates for a fluid with a tilted surface knowing 
only the boundary function and the tilt angle. The first step is to compute 
the height (axial value) of the point on the low end of the surface as shown 
in the following illustration: 



( 1 ) 

( 2 ) 

The value Zg£^ is obtained by solving this linear equation simultaneously 
with r = f(Z). Once these values are found, coordinates are generated at 
each elevation, circumferential and radial station by interpolation and solu- 
tion of Eq. (2). At each step, GRID Bulk Data is generated. A subsequent 
call to SEL0UT generates the element connection cards. 


The equation of the surface line is defined by; 



2 = mr + Zj^ - mf(Zj^) 
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7.5 SHLCAP 


7.5.1 Entry Point : SHLCAP 

7.5.2 Purpose 

To generate both grid points and element connections for the shell 
cap model. 

7.5.3 Calling Sequence 
CALL SHLCAP 

C0MM0N/C0NTRL/ C0MM0N/ELTAB/ 

C0MM0N/MESH/ C0MM0N/PSDC/ 
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7.6 S0LCAP 

7.6.1 Entry Point ; S0LCAP 

7.6.2 Purpose 

To generate both grid points and element connections for the solid cap 
model. In addition, to generate PL0TELs for the fluid free surface, if 
applicable. 


7.6.3 Calling Sequence 

CALL S0LCaP 

C0MM0N/MLINE/ 

C0MM0N/BTAB/ 

C0MM0N/SHP/ 


C0MM0N/C0NTRL/ 

C0MM0N/MESH/ 

C0MM0N/ELEM/ 



7.7 STRNGR 


7.7.1 Entry Point : STWiGR 

7 . 7.2 Purpose 

To generate stringer connections (CBAR or CR0O) defined with shell models. 

7.7.3 Calling Sequence 

CALL STRNGR 

C0MM0N/MESH/ 

C0MM0N/MLINE/ 
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7.8 FLSTRl 


7.8.1 Entry Point ; FLSTRl 

7.8.2 Purpose 

To determine the fluid/structure interfaces and generate CFLSTR Bulk Data 
cards. 

7.8.3 Calling Sequence 
CALL FLSTRl (IC0DE) 

where : 

IC0DE - specifies whether the shell elements are on the SAVE file (IU3) 
or the scratch file (IU4) 

7.8.4 Method 

The structure boundary ranges (TMIN, TMAX, 2MIN, ZMAX; see Section 2.1.3) are 
read from IU3 or IU4 into core. The fluid element ranges are then read from 
the SAVE file one at a time and compared with the shell ranges. All overlapping 
shell element ID's are saved. When the process is completed for a given fluid 
element, a CFLSTR Bulk Data card (or cards) is generated and the next fluid 
element is read. This process continues until all fluid elements have been 
tested. 
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8. SYSTEM CONSIDERATIONS 


MESHGEN was written for, and Implemented on, an IBM 360/165. It is currently 
being implemented on the UNIVAC 1108. With the overlay structure defined in 
Figure 3 , MESHGEN will execute in 180K decimal bytes (45K words) . This 
number can be affected by local modifications to system routines, FORTRAN 10 
packages, or SC4020 software. 

A list of subroutines that can cause implementation problems when installing 
on a new system are summarized. 

8.1 OVERLAY FOR MESHGEN 


ROOT 

SEGMENT 


OVERLAY B 


CCARD 

STR0UT 

WED0UT 

TET0UT 

HEX0UT 

PL0TEL 

FRE0UT 

GETLIN 



PL0T2 

PLTNFR 

PLTRD 

TPLAB 

TPLIN 

TYPLl 

TYPL2 

TYPL3 

TYPL4 


XRCARD 

LEXER 

LEVELl 

LEVEL2 

LEVEL3 

LEVEL4 

LEVELS 


MAIN 

GRD6EN 

BUG 

ZER0UT 

MESHBD 

FUNC 

MAPFNS 

PAGE 

K L 

INT2AL 

DATAEX 

PLTHED 

INIT 

GETT0K 

DIVLIN 

ENC0DE 

INTERP 

GETCRD 

PIFl 

PNCHER 

DEC0DE 

MESAGE 

STRNGR 

WRTHED 

GRD0UT 

PREPAS 

GETVAL 

APRXEQ 



FLSTRl 

CHECKR 


FIGURE 3. OVERLAY FOR MESHGEN 
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8.2 MACHINE OR FACILITY DEPENDENT ROUTINES 


8.2.1 MAPFNS 

This routine is written In assembly language and must be written for each 
system before Implementation. 

8.2.2 PREPAS 

Used to convert EBCDIC to BCD before processing begins. For systems utilizing 
BCD only, the routine should be rewritten to simply copy the Input. 

8.2.3 MESHBD 

The machine constants defined In this BLOCK DATA should be changed to reflect 
the system architecture. 

8.2.4 PL0TS 

The plot driver routine, PL0T2, as well as all plot interfaces should be 
checked to conform to facility-dependent conversion. 
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